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The history of the Universe after its first second is now tested by high quality observations of light element 
abundances and temperature anisotropies of the cosmic microwave background. The epoch of the first 
second itself has not been tested directly yet; however, it is constrained by experiments at particle and heavy 
ion accelerators. Here I attempt to describe the epoch between the electroweak transition and the primordial 
nucleosynthesis. 

The most dramatic event in that era is the quark-hadron transition at 10 /is. Quarks and gluons condense 
to form a gas of nucleons and light mesons, the latter decay subsequently. At the end of the first second, neu- 
trinos and neutrons decouple from the radiation fluid. The quark-hadron transition and dissipative processes 
during the first second prepare the initial conditions for the synthesis of the first nuclei. 

As for the cold dark matter (CDM), WIMPs (weakly interacting massive particles) — the most popular 
candidates for the CDM — decouple from the presently known forms of matter, chemically (freeze-out) 
at 10 ns and kinetically at 1 ms. The chemical decoupling fixes their present abundances and dissipative 
processes during and after thermal decoupling set the scale for the very first WIMP clouds. 
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1 Introduction 

This article is an attempt to summarize our present understanding of the early Universe in the epoch be- 
tween the electroweak (EW) transition and the onset of primordial nucleosynthesis. During that time the 
temperature drops from Tew ~ 100-200 GeV at the EW transition to 1 MeV, when all weak interaction 
rates fall below the Hubble rate, which marks the beginning of the primordial nucleosynthesis epoch. At 
1 MeV, the Hubble age, tn = 1/H, is 1 s; in that sense this is a review about the first second of the Uni- 
verse. The physics of the first second is well tested at high energy colliders, as far as the particle content 
up to masses of ^ 100 GeV is concerned. The equation of state of the radiation fluid that dominates the 
Universe during the first second is not known up to the highest temperatures, but it is tested cuiTently by 
heavy ion experiments up to the scale of quantum chromodynamics (QCD) at 100-200 MeV. However, I 
would like to stress that, since details of the QCD transition and the nature of cold dark matter (CDM) are 
still unknown, we have to rely on yet untested models. 

The Universe was radiation-dominated during the first second. Today, this conclusion can be based on 
several lines of argument. The most direct evidence is the observed rise of the temperature T of the cosmic 
microwave background (CMB) radiation with increasing redshift z, i.e. T{z) = To(l + z) fl\. Together 
with the observed Planck spectrum of the CMB (with Tq = 2.725 ± 0.001 K (2|'), we can conclude that 
the radiation energy density is given by the Stefan-Boltzmann law, 

e=^5.(T)T4(x(l + z)^ (1.1) 

which, at large enough z, dominates the energy density of matter Cm = nm {l + z)^. The function ge{T) 
counts the effective number of relativistic helicity degrees of freedom at a given photon temperature T 

' In the following we will set Boltzmann's constant k = 1, and thus measure temperature in units of eV, i.e. 1 K = 8.617x 
10^^ eV. Unless stated otherwise, we set c = ?i = 1. 
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(fermionic degrees of freedom are suppressed by a factor 7/8 with respect to bosonic degrees of freedom). 
For times after e+e~ annihilation (at T^+g- ^ me/3 « 170 keV) one finds 

i = 2.™xl0-(^)(f).-l±i, (1.2) 

where measures the mass density of matter^; the redshift of matter-radiation equality is given by 
1 + Zcq = 3612(cjm/0.15). Before e+e^ annihilation we have 

i ^ — «»(a(Si)(Tl^) 

= 0.830 fgPU^I; (1.4) 



10.75 J V 1 + Zcq 



1 + . = 5.966 X 10^ . (1.5) 

Thus, during the first second, radiation is dominating matter and all other possible components of the 
Universe (such as curvature, which scales with a^^, or a cosmological constant, which does not change 
during the cosmic evolution). 

It is very useful to introduce the fundamental cosmological scale, associated with a given temperature 
of the Universe. The Friedmann equation links the expansion (Hubble) rate H to the mass density of the 
Universe, 

H^='-^P. (1.6) 

where p = tjf?. This equation is obtained under the assumption that general relativity is the appropriate 
description of gravity and that the space-time is isotropic and homogeneous, with vanishing spatial curva- 
ture and cosmological constant. For dynamical questions, the Hubble time"', = 1/ff, is the typical time 
interval in which any of the thermodynamic variables, the curvature and the expansion of the Universe, can 
change significantly: 

/10.75\^^^ /lMeV\^ 

The time scale of the QCD transition is 10 /is, that of the EW transition is 10 ps. It is straightforward to 
convert these time scales into distances; the Hubble distance i?H = cin at 1(160, 10^) MeV is 4 x 10^ km 
(10 km, 10 mm). These distances are the physical distances at the given temperatures and they grow with 
the expansion. 

Let us now consider the effective number of relativistic helicity degrees of freedom as a function of the 
temperature, see fig. II. II The two full lines show the effective degrees of freedom of the energy density, 

and of the entropy density (upper line) for the particle content of the standard model. In this plot 
the measured particle masses are taken from the Particle Data Group |5|. With increasing temperature the 
following important events happen: at temperatures below T^+g- ~ me/3, present values ge(T'o) = 3.363 
and gs{To) — 3.909 are taken. The difference between g^ and gs, the effective number of relativistic 



T 



Ljjn measures the mass of non-relativistic matter in the Universe today, i.e. baryons and cold dark matter It is defined as 
oJm = Pm/[3(i?o//i)^/87rG], where pm denotes the mass density and the expansion (Hubble) rate is given by Hq = 
lOO/ikm/s/Mpc. Recent data from WMAP, combined with data from other CMB experiments and the 2dF galaxy redshift 
survey, provide = 0.136 ± 0.009 and h = 0.71 ± 0.04, within a fit to a ACDM model with running spectral index (^J- 
In a radiation-dominated Universe the cosmic time is proportional to the Hubble time, t = However, in this work 

we refer to the Hubble age, since only lower limits can be given for the age of the Universe; e.g. a past infinite epoch of 
cosmological inflation (avoiding any singularity) is possible before the radiation-dominated epoch 13 • 
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MSSM 




Fig. 1.1 The effective number of degrees of freedom (?£ (T) = e(r)/(7r^/30r'*). The fullline is the prediction of the 
standard model of particle physics, the dashed line shows a minimal supersymmetric extension of the standard model 
(SPSla)L6J. 



helicity degrees of freedom that contribute to the entropy density, is due to the difference of the photon 
and neutrino temperatures T^/T^ = (4/ll)^/'^ after e+e^ annihilation. At temperatures above 1 MeV, 
electrons, photons and neutrinos have the same temperature. As the temperature increases, particles with 
mass m become relativistic at approximately m/S. The rise of g starting at around 30 MeV is mainly due 
to muons and pions, but also heavier hadrons can be excited. In figure [TTTI all hadrons up to the uj meson 
(m„ = 783 MeV) have been considered (finite volume effects are ignored). At the temperature of the 
QCD transition, here Tqcd — 160 MeV, the number of degrees of freedom changes very rapidly, since 
quarks and gluons are coloured. As discussed in more detail below, the order of the QCD transition is 
still unknown. Thus the details of figure FTTTl in the vicinity of the QCD transition are just an indication of 
what could happen. The strange mass is assumed to be 120 MeV here. At still higher temperatures again 
heavier particles are excited, but within the standard model of particle physics nothing spectacular happens. 
Especially the EW transition is only a tiny effect in figure FTTTI (Here, we assumed a top mass of 174 GeV 
and a Higgs mass of 120 GeV.) This situation changes if the minimal supersymmetric extension of the 
standard model (MSSM) is considered. The effective degrees of freedom are more than doubled, since all 
supersymmetric partners can eventually be excited"^. For the purpose of figure FTTTl the mass spectrum of the 
"Snowmass Points and Slopes" model la |6| has been used. Note that this model has light superpartners, 
compared with the other models. The lightest supersymmetric particle in the SPSla model is the lightest 
neutralino with a mass of 95 GeV. 

Figure [LTI clearlv highlights three potentially very interesting epochs: (i) the EW transition at 100- 
200 GeV, (ii) the QCD transition at 150-180 MeV, and (iii) the e+e" annihilation at 170 keV. At these 
events the number of relativistic helicity degrees of freedom changes dramatically, and so the entropy 
density of the radiation fluid. Figure [1~T] is therefore very useful to tell us which events in the radiation 
fluid might be of interest from the thermodynamic point of view. 

e^e'^ annihilation is treated in detail in many cosmology textbooks and we refer the reader to (TJ. 
The three neutrino flavours are decoupled chemically and kinetically from the plasma at temperatures 
below 1 MeV; thus the entropy of the relativistic electrons is transfeiTed to the photon entropy, but not 
to the neutrino entropy when electrons and positrons annihilate. This leads to an increase of the photon 
temperature relative to the neutrino temperature by TyjT^ — (4/11)^/"^. 



The standard model has g = 106.75, the MSSM has g = 228.75, when the temperature is larger than all particle masses. 
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Fig. 1.2 The Hubble rate H and the typical interaction rates of weak (r„) and electric (Fe) processes that involve 
relativistic particles, as well as the typical rate of a weak annihilation rate r„.ann for a particle mass of 100 GeV. At 
~ 1 s, the weak interaction rate falls below the expansion rate (chemical and kinetic decoupling of neutrinos, kinetic 
decoupling of neutralinos), at temperatures of the order of 1-10 GeV neutralinos freeze-out. The electric interaction 
rate stays well above the Hubble rate up to the epoch of photon decoupling, well after the epoch we are interested 
in here. The dashed line indicates the rate 1/s. (Note that for the purpose of this sketch the weak interactions are 
approximated by the Fermi theory, which breaks down at T '-^ mw/S. We also ignored the running of the coupling 
constants in this figure.) 

For the EW transition we still do not know the particle content of the Universe at the relevant tempera- 
ture. We will not study its cosmological implications here. During this transition, according to the standard 
model of particle physics, all particles except the Higgs acquire their mass by the mechanism of sponta- 
neous symmetry breaking. For a Higgs mass above 72 GeV (the current mass limit is 114.3 GeV at 95% 
CL |3|), it has been shown that the transition is a crossover, and since the change in relativistic degrees of 
freedom is tiny, this is a very boring event from the thermodynamical perspective. However, extensions of 
the standard model offer the possibility that the transition is of first order, which implies that electroweak 
baryogenesis might be possible in such a case (although the region of parameter space in which this can 
happen is small). For recent results on the order and temperature of the electroweak transition, see 1 8 1. A 
cosmologically most relevant consequence of the epoch before the electroweak transition is that sphaleron 
processes allow a violation of baryon number B and of lepton number L, in a way that A(i? — L) = 
[9J. At high temperatures T > Tew; these processes take place at a rate per unit volume of the order 
g^^ \n.{l /g)T'^ (g denotes the SU(2) coupling) 1101 . which leads to a rapid equilibration between baryon 
and lepton number This has the important consequence that (within the standard model) the sphaleron 
processes predict L = —(51/28)5. The numerical coefficient depends on the particle content before the 
electroweak transition, e.g. for two complex Higgs doublets L = —(15/8)5 II II . In any case the lepton 
asymmetry has to be of the same order of magnitude as the baryon asymmetry. This is important for the 
question whether chemical potentials can play a role in the early Universe. We show below (section l3r4l 
that they are completely irrelevant for the present purpose. 

Thus, from the thermodynamic point of view, the QCD epoch is the most interesting event during the 
first second. However, not all interesting physics is covered by thermodynamics. Especially phenomena 
that affect the baryon content and cold dark matter are not only governed by the bulk thermodynamics. We 
therefore take a closer look at the interaction rates. 

A complementary picture of the early Universe arises from a comparison of the most relevant interaction 
rates and the expansion rate, sketched in figure n~2l as a function of temperature. Strong, electric, and weak 
interactions keep all relativistic particles in kinetic and chemical equilibrium down to temperatures of 
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~ 1 MeV. At that point neutrinos and neutrons decouple chemically and kinetically from the rest of the 
radiation fluid. This has several important implications: (i) During the process of neutrino decoupling, 
density inhomogeneities of the radiation fluid on scales below 1/30 of the Hubble distance at 1 MeV are 
washed out by collisional damping. Irrespectively of the initial conditions of the Universe, this guarantees 
that entropy is distributed homogeneously within large patches, and thus justifies neglecting temperature 
fluctuations on small scales during big bang nucleosynthesis (BBN). (ii) The weak equilibrium between 
neutrons and protons freezes out as the neutrinos decouple. The neutron-to-proton ratio changes only 
because of neutron decay after weak decoupling, (iii) Neutrons no longer scatter and can travel "large" 
distances. They can smear out local fluctuations of the neutron-to-baryon ratio. 

Because of neutrino oscillations, the individual neutrino flavour is not conserved, but a flavour equilib- 
rium is established at a time scale of the order of the oscillation time, as soon as the radiation fluid is dilute 
enough to allow a sufficient mean free path for the neutrinos. For a scenario based on the atmospheric 
and solar neutrino data, favouring large mixing angles, it was shown by Dolgov et al. 1121 that flavour 
equilibrium is established at about T ^ 3 MeV. 

Weakly interacting massive particles (WIMPs) are excellent candidates for the CDM. A prominent 
example is the lightest supersymmetric particle (LSP), most likely the lightest neutralino. WIMPs are 
exponentially suppressed as the temperature drops below their mass. However, since they interact only 
weakly, the annihilation rate cannot keep up with the Hubble rate and they drop out of chemical equilibrium 
(freeze-out). This happens at typically T ^ m/20. In figure [O] m = 100 GeV. Nevertheless, they are 
kept in kinetic equilibrium down to temperatures of 10-100 MeV |T3l by elastic scattering. During and 
after their kinetic decoupling, collisional damping and free streaming wash out density inhomogeneities in 
CDM on very small scales. This effect might be of relevance for structure formation and the search for 
dark matter. 

The notation used in this review is summarized in two tables at the end of the paper. 

2 The cosmic QCD transition: an overview 

One of the most spectacular epochs in the early Universe is the QCD epoch. QCD describes the strong 
interactions between quarks and gluons and is well tested in the perturbative regime, i.e. at high energies 
and momenta. At low energies, quarks and gluons are confined into hadrons. The scale of QCD is Aqcd = 
0(100 MeV). At temperatures T Aqcd, there is a transition between a quark-gluon plasma (QGP) and 
a hadron gas (HG) 114|. These theoretical expectations are in agreement with findings from the CERN 
heavy ion programme II15I and the ongoing studies at the Relativistic Heavy Ion Collider (RHIC) 1161 at 
BNL. 

The first studies of the cosmological QCD transition date back to the early 80's 1171 . It was then reaUzed 
that a first-order QCD transition proceeds via bubble nucleation I18II19I with small supercooling I20II21I . 
The bubbles were shown to grow most probably by weak deflagration II22II23I . Most of the following work 
has been based on this scenario of a first-order QCD transition with homogeneous bubble nucleation and 
bubble growth by weak deflagration. A detailed discussion of this scenario has been given by Kajantie and 
Kurki-Suonio |24|. 

Based on this scenario a separation of cosmic phases has been suggested by Witten 1201 . which may give 
rise to large inhomogeneities in the distribution of baryons in the Universe. Originally it was hoped that 
the mean bubble nucleation distance is a reasonable fraction of the size of the Universe, the Hubble scale. 
From recent lattice QCD results for the latent heat and surface tension in the quenched approximation (no 
dynamical quark flavours) 125 1 , a tiny mean bubble nucleation distance follows 1261 1271 l28l . 

The mean bubble nucleation distance might be larger for other bubble nucleation scenarios. Nucleation 
at impurities has been investigated by Christiansen and Madsen |28|. Another possibility is that the nu- 
cleation of bubbles is determined by pre-existing inhomogeneities of temperature (e.g. generated during 
inflation) 1 29 1, or by pre-existing baryon number inhomogeneities 1301 . 
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Apart from the question 'What happened at the cosmological QCD transition?', the cosmological QCD 
transition is of interest for two more reasons: firstly, the QCD epoch sets the initial conditions for one of 
the pillars of standard cosmology — big bang nucleosynthesis. The QCD transition and baryon transport 
thereafter determine the distribution of protons and neutrons at the beginning of nucleosynthesis. It was 
first recognized by Applegate and Hogan |21| that an inhomogeneous distribution of baryons due to a 
first-order QCD transition could change the primordial abundances of the light elements. This issue has 
attracted a lot of interest 1311 l32l 1331 1341 because inhomogeneous BBN might allow a higher baryon-to- 
photon ratio than homogeneous BBN [351. However, today this possibility is ruled out by observations of 
CMB temperature anisotropics 1,36,.. 37. . 3J. which confirm independently the value of cuh (defined as Wm, 
but for the baryon mass density) that is extracted from the theory of homogeneous BBN and the observed 
light element abundances^. However, with the increasing precision of the determination of abundances, we 
might find ourselves in a situation in which some tension arises between the values of ujt, as extracted from 
the individual measurements of "^He, D, and ^Li; see I39II40I for an indication that there might be a conflict 
between the observed D abundance, which is consistent with the CMB measurements, and ^He and ^Li 
measurements. Besides systematic errors of the abundance determination, this might be a signature for an 
inhomogeneous BBN. Moreover, an inhomogeneous BBN might allow the primordial production of heavy 
elements such as ^^C |41 1. 

Secondly, the QCD transition might generate relics that might be observable today. Most of these relics 
are only formed if the QCD transition is first-order and if the latent heat is much larger than that found from 
quenched QCD in lattice calculations. Strange quark nuggets as dark matter and gravitational waves from 
colliding bubbles have been suggested by Witten |20|. Hogan [42 1 and more recently Cheng and Olinto 
1431 argued that magnetic fields might be generated. Today, these relics appear to be unlikely, because the 
typical scale for bubble nucleation is very small. 

Recently, Zhitnitsky [44| suggested a new CDM candidate: QCD balls. If axions existed and if the 
reheating scale after inflation lied above the Peccei-Quinn scale, collapsing axion domain walls could trap 
a large number of quarks. At some point the collapse would be stopped by the Fermi pressure of the 
quarks, which would settle in a colour superconducting phase |45|. This process takes place during the 
QCD transition, but does not require a first-order transition, contrary to the idea of strange quark nuggets. 
Brandenberger, Halperin, and Zhitnitsky | 46 | speculated that even a separation of baryons and antibaryons 
due to the non-perturbative QCD vacuum might be possible during the QCD epoch. This would effectively 
give a new baryogenesis scenario and a new candidate for dark matter. It seems to me that much more 
work has to be done to explore these exciting ideas. According to the mentioned suggestions, in the best 
of all QCD worlds, it might be possible to explain baryogenesis and the nature of CDM! 

It was speculated for various reasons that black holes could form during the QCD transition I47ll48ll49l . 
From the present-day perspective these scenarios seem to be highly unlikely. 

The QCD transition itself might also lead to the formation of small CDM clumps I50ll51ll52l . The speed 
of sound vanishes during a first-order QCD transition and thus the restoring forces vanish. This leads to 
large amplifications of primordial density fluctuations in the radiation fluid and in cold dark matter. This 
mechanism can work for axions or primordial black holes, since they are kinetically decoupled at the QCD 
transition, but not for WIMPs, as they belong to the radiation fluid during the transition. 

Independently from the order of the transition, a primordial background of gravitational waves is mod- 
ified by the QCD transition, as shown in |53|. 

Reviews on the cosmological QCD transition may be found in Refs. I54ll55ll56ll57ll58l . 

2.1 Scales 

The QCD transition is expected to take place at ^ Aqcd- We would Uke to measure from heavy ion 
colUsions, but it turns out that this is not a simple task. A temperature estimate can be obtained from the 



^ ojb = 0.022 ± 0.001 from a seven-parameter fit to CMB and large scale structure data (3], whereas o^b = 0.0205 ± 0.0018 
from the most recent determination of the primordial deuterium abundance 1381 . 
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study of measured hadron abundances in reletivistic heavy ion collisions. After some short initial phase 
one expects that a quark-gluon plasma (QGP) is formed, which expands and eventually makes a thermal 
confinement transition. At some point the hadron gas is dilute enough such that the abundance of hadron 
species is fixed. At RHIC experiments, this tempertature of the chemical decoupling of hadrons has been 
estimated to be T = 174± 7 MeV at a baryon chemical potential /ib = 46± 5 MeV (statistical errors only) 
L59J . This reasoning suggest that the QCD transition temperature should lie above the estimated hadron 
freeze-out temperature and thus T^, > 160 MeV (95% C.L.). Although the baryon chemical potential 
at RHIC is small with respect to previous heavy ion experiments, it is huge from the point of view of 
cosmology. However, lattice QCD predicts a suprisingly small curvature of the transition temperature as 
a function of the baryon chemical potential at /ie = (for 2 + 1 quark flavours, see 1601 1. Thus, the 
lower limit on the transition temperature suggested by RHIC should be applicable in cosmology. Another 
difference between heavy ion collisons and cosmology should be stressed: while the time scale of the 
cosmological QCD transition is 10^^ s, it is lO^^'^ s in the laboratory. It is therefore necessary to check 
the theoretical and experimental estimates of by lattice QCD 'experiments'. 

From recent lattice QCD calculations for quenched QCD (no dynamical quarks) the transition tempera- 
ture is « 271 ± 2 MeV |6l]|62l|6a. For two-flavour QCD « 171 ± 4 MeV |6l|65l|66||6a, whereas 
for three-flavour QCD « 154 ± 8 MeV |64, 63], almost independent of the quark mass. For the most 
interesting case of two light quark flavours (up and down) and one massive strange quark, no values for 
the transition temperature have been obtained so far. In the following I pick a transition temperature of 

= 160 MeV, i.e. I implicitly assume that the physical situation resembles more closely the three-flavour 
case than the two-flavour situation. Recent reviews of thermal QCD simulations can be found in Ii63il . 



2.1.1 The Hubble sc ale 

It is a good approximation for our purpose to treat all particles with m <^ 3T as though they were mass- 
less; all other particles are neglected in the total energy density e. Above the QCD transition ^quarks = 
(7/8)12iV/ (number of quark flavours) and (jgiuons — 8, below the QCD transition ^pions = 3. At the QCD 
epoch there are the photons, three flavours of neutrinos, and electrons and muons. Counting all up we find 
for the QCD epoch 

5(T > T*) = 51.25 (61.75) , (2.1) 

without (with) strange quarks, and 

g(T < T*) = 17.25 (21.25) , (2.2) 

without (with) kaons. 

At the QCD transition the Hubble radius is about 10 km; 



61.75y /160MeV\ 



2 



Rh - [ j ) 7.2 km , (2.3) 



before the transition, and 



i?H ~ ( ^ ) ( ^""J" ' ) 12km (2.4) 



21.25y 1^160 MeVy 



after it. 

Today this corresponds to scales of 1 pc or 3 light-years. The Hubble time at the QCD transition, 
^ 10^^ s, is extremely long in comparison with the relaxation time scale of the strong interactions, 
which is about 1 fm/c ~ 10"^'^ s. Thus, the transition is very close to an equilibrium process. 
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a- 



a+ 



Fig. 2.1 Qualitative behaviour of the temperature T as a function of the scale factor a during a first-order QCD transi- 
tion with small supercooling. Above the critical temperature T* the Universe (filled with a quark-gluon plasma) cools 
down thanks to its expansion (a < a_). After a tiny period of supercooling (in the figure the amount of supercooling 
and its duration are exaggerated) bubbles of the new phase (hadron gas) nucleate. Almost in an instant these bubbles 
release enough energy (latent heat) to reheat the Universe to the critical temperature. During the rest of the transition 
both phases coexist in pressure and temperature equilibrium (a_ < a < a+). Therefore the temperature is constant. 
After reheating the bubble growth is governed by the expansion of the Universe. The transition is completed when no 
quark-gluon plasma is left over. For a > a+ the temperature decreases again due to the expansion of the Universe. 



The mass inside a Hubble volume is ~ IM, 



0- 



61.75~\ = i^ieOMeVy 



2. 5 Mr, 







(2.5) 



This mass is redshifted cx (1 + z) as the Universe expands, because it is made up of radiation. An invariant 
mass is the mass of cold dark matter in a comoving volume, Afcdm = (1 + Zcq) /(I + z)M{z). At the QCD 
transition (1 + Zeq)/{1 + ~ 10~®, and thus 







(2.6) 



Another figure of interest is the baryon number inside the Hubble volume at the QCD transition, Bh = 
(47r/3)i?^nB. The baryon number density ub follows from the ratio of baryons to photons at BBN, 
?7 = (ne /?^7)bbn, and the conservation of baryon number and entropy, ub/s = constant: 



S / BBN 



(2.7) 



Finally, the baryon number inside a Hubble volume reads: 

3 



^ ^ ^ 61.75 \ / 160MeV \ 



9 



V 



/ V 5.6x10 



2.1 X 10' 



48 



(2.8) 



at the beginning of the transition and about twice that value at the end . From standard BBN and the most 



recent deuterium measurements, rj = (5.6 ± 0.5) x 10 13^ 



This formula is correct if no black holes are formed during the QCD transition and if the quark nuggets that might have formed 
evaporate before the BBN epoch. 
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2.1.2 The bubble scale 

The typical duration of a first-order QCD transition is O.ltn- If the cosmological QCD transition is of first 
order, it proceeds via bubble nucleation. From the small values of surface tension and latent heat found in 
lattice QCD calculations |25 1, the amount of supercooling is found to be small [27 1. The temperature as 
a function of the scale factor is shown in fig. l2.1l for the small supercooling scenario. Hadronic bubbles 
nucleate during a short period of supercooling, At^c ^ 10^'^<h- The typical bubble nucleation distance is 

dnuc - 1 cm - 10-^i?H (2.9) 

for homogeneous nucleation |28 1. The hadronic bubbles grow very fast, within 10~*'tH, until the released 
latent heat has reheated the Universe to T^,. By that time, just a little fraction of volume has gone through 
the transition. For the remaining 99% of the transition, the HG and the QGP coexist at the pressure 
Phg(7*) = PQGp(T'it)- During this time the hadronic bubbles grow slowly and the released latent heat 
keeps the temperature constant until the transition is completed. The energy density decreases continuously 
from eQGp(T*) at the beginning of the transition to eHG(2*) when the transition is completed. 



2.2 Order of the thermal QCD transition 

The order of the QCD phase transition and (for a first-order transition) the magnitude of the latent heat are 
still a subject of debate. Below I first introduce the bag model, because it gives a simple parametrization 
for the pressure p{T) and it is useful to introduce some notation. Then I sum up the knowledge obtained 
from lattice QCD (order of the transition, latent heat, surface tension). 



2.2. 1 The bag model 

The MIT bag model 1167 1 represents the short-distance dynamics by an collisionless gas of massless quarks 
and gluons [in the following I call that a Stefan-Boltzmann (SB) gas for brevity] and the long-distance 
confinement effects by a constant negative contribution to the pressure, the bag constant B, 

PQGP(r)=pSBp(r)-B, (2.10) 

withpQQp(T) = (7r^/90)5QGpT''*. Here I take gqcp = 47.5. The low-temperature phase is a hadron gas. 
It may be modelled as an SB gas of massless pions and kaons (since ttik ~ 3T*,), ^hg — 7, 

Phg = -^gacT^- (2.11) 
At the phase transition the pressures of the quark-gluon phase and the hadron phase are in equilibrium, 
PQGp(r*) = Pbg{T*) = . (2.12) 
This condition gives, together with ( I2.10> and ( I2.11> . the relation between and B: 

S=^(.9QGP-5HG)T^ (2.13) 

Let me take T^, from lattice QCD calculations with two and three flavours, which indicate = 150 to 
180 MeV |63 1. This corresponds to a range of bag constants S^/^ = 218 to 261 MeV. This range overlaps 
well with fits to the light-hadron masses, which yield S^/** = 145 to 245 MeV (a compilation of various 
bag model light-hadron fits can be found in |68 1). 

The QCD transition is first order in the bag model, because the entropy density s = dp/dT makes a 
jump. The latent heat per unit volume 



I EE T^As 



(2.14) 
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measures the amount of 'internal' energy that is released during the phase transition. In the bag model the 
latent heat is 

I = -^AgT^ = 45, (2.15) 

where Ag = ggcp - 5hg- 

Besides the latent heat, the surface tension is the crucial parameter for the nucleation of bubbles (see 
section|5}- The surface tension 




(2.16) 



is the work dW that has to be done per area dA to change the phase interface at fixed volume. The absence 
of surface excitations in hadronic spectra suggests that a^^^ ^ B^^^ 1681 in the bag model. This implies 
fj ^ 3.irf for three quark flavours. A self-consistent calculation of the surface tension within the MIT 
bag model shows that the surface tension vanishes for massless quarks and gluons if no interactions besides 
the bag constant (i.e. as = 0) are taken into account II68I . This can be cured by introducing short-range 
interactions f681 or by including the strange quark with mass 1691 . nis — 80 to 150 MeV |5 1. 



2.2.2 Lattice QCD results 
In quenched QCD the phase transition is of first order 1611 . The latent heat was determined to be Z w I 'iT'^ 



1251 . It is useful to take the ratio Rl of the latent heat to the value As^^, where As^^ is the difference 
in entropy between two SB gases: 

I _\ 1 bag model 
(r*As)SB \ 0.2 quenched lattice QCD ' ^ ' 

For two light quarks it is likely that the transition is a crossover 1701 1711 . This is in agreement with 
theoretical considerations [72], which predict a second-order phase transition in the massless quark limit. 

For three flavours close to the chiral limit, the phase transition is again of first order and some simula- 
tions suggest that this holds true for the physical case 1731 . The latter result was obtained using the Wilson 
quark action, whereas results with staggered quarks f74' '641 indicate a crossover for the physical quark 
masses. For four quark flavours the transition is first order L75,l . For a detailed reviews on these issues, see 
|63. 

Since the latent heat for lattice QCD is known for quenched QCD only, I decided to use the latent heat 
ratio Rl = 0.2 from quenched QCD as an indication for the physical case. 

For the surface tension rather small values are found from lattice QCD. It is reported in |25 1 that a « 
O.OIST^ for quenched lattice QCD. There are no values for unquenched QCD available yet. However, an 
upper bound was obtained for the case of four-flavour lattice QCD in Ref. li76J . i.e. a < O.IT^. 



2.3 Effects from a first-order QCD transition 

Let me now briefly summarize the effects that have been suggested to emerge from the cosmological QCD 
transition. There are two kinds of effects: the effects that have been found in the mid 80s and early 90s stem 
from the bubble scale and they thus affect scales A < c?nuc- The formation of quark nuggets, the generation 
of isothermal baryon fluctuations, the generation of magnetic fields and gravitational waves belong to the 
effects from the bubble scale. 

In recent years it was found that there is another class of possible consequences from the QCD transition, 
which are connected to the Hubble scale and therefore affect scales A < i?H. Among these effects are 
the amplification of inhomogeneities and later formation of cold dark matter clumps, the modification 
of primordial gravitational waves, and the enhanced probability of black hole formation during the QCD 
transition. 
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2.3.1 Quark nuggets/Strangelets 

In the mid 80s interest in the cosmological QCD transition arose, because it was realized that a strong 
first-order QCD phase transition could lead to important observable signatures. Most of the interest was 
based on a separation of cosmic phase as suggested by Witten |20|. If the cosmological QCD transition is 
first-order, bubbles of hadron gas are nucleated and grow until they merge and fill up the whole Universe. 
Towards the end of the coexistence of the QGP and HG phases, shrinking quark droplets remain. These 
droplets are expected to be baryon-enriched with respect to the hadron phase since in equilibrium, baryons 
are suppressed by the Boltzmann factor in the HG, and baryon diffusion across the phase boundary might 
be inefficient. 

In 1971 Bodmer |77| suggested the possibility that strange quark matter might be the ground state of 
bulk matter, instead of '^^Fe. Later Witten 1 20 1 discovered this idea again. Strange quark matter was further 
studied by Farhi and Jaffe |68 1. The idea of strange quark matter is based on the observation that the Pauli 
principle allows more quarks to be packed into a fixed volume in phase space if three instead of two flavours 
are available. Thus the energy per baryon would be lower in strange quark matter than in nuclei. However, 
the strange quark is heavy compared with up and down quarks, and this mass counteracts the advantage 
from the Pauli principle. No strange quark matter has been found experimentally so far 1781 . The issue of 
stability of strange quark matter has not been settled yet; for a recent review see 1791 . 

Witten 1201 pointed out that a separation of phases during the coexistence of the hadronic and the 
quark phase could gather a large number of baryons in strange quark nuggets |20|. These quark nuggets 
could contribute to the dark matter today |20| or affect BEN |80|. At the end of the transition the baryon 
number in the quark droplets could exceed the baryon number in the hadron phase by several orders of 
magnitude, could be close to nuclear density II81I . However, it was realized that the quark nuggets, 

while cooling, lose baryons. The quark nuggets evaporate as long as the temperature is above ^ 50 MeV 
1 82 1. Quark nuggets may survive this evaporation if they contain much more than ^ 10'''* baryons initially 
1 83 1. This number should be compared with the number of baryons inside a Hubble volume at the QCD 
transition, which is 10''^ (see section |2.1.1> . Thus, the mean bubble nucleation distance should be > 
3 X 10~^i?H ^ 300 m so as to collect enough baryons. This seems impossible from recent lattice QCD 
calculations of latent heat and surface tension l25l . 

In I81II83I a chromoelectric flux-tube model was used to estimate the penetration rate of baryons through 
the interface. A quark that tries to penetrate the interface creates a flux tube, which most probably breaks 
up into a quark-antiquark pair. By this mechanism, mesons evaporate easily. On the other hand, baryons 
are rarely formed, because a diquark-antidiquark pair has to be produced in the break up of the flux 
tube. However, one could think of mechanisms that would increase the evaporation rate of baryons. If a 
significant fraction of diquarks was formed in the quark phase, these diquarks could penetrate the interface 
by creating a flux tube, which eventually breaks, creating a quark-antiquark pair. Then the quark would 
evaporate together with the diquark and form a baryon, whereas the antidiquark would remain in the quark 
phase. Such a mechanism would increase the evaporation rate. Thus, independently from the existence 
(stability) of strange quark matter it seems highly unlikely that strange quark nuggets could survive after 
the cosmological QCD transition below temperatures of ^ 50 MeV. 

However, see Bhattacharyya et al. |84| for a different point of view. They conclude that stable quark 
nuggets could be formed if the QCD transition is a strong first-order transition (they use the bag model) and 
if the critical temperature is around 100 MeV, instead of 150-180 MeV as indicated by lattice QCD. They 
speculate that these quark nuggets could account for all the CDM. In a recent work [85], it was suggested 
that these primordial quark nuggets might clump by gravitational attraction and eventually form bound 
objects of O.5M0, which would explain the gravitational microlensing events that have been observed 
towards the Large Magellanic Cloud 1861 . 
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2.3.2 Inhomogeneous nucleosynthesis 

Applegate and Hogan 1211 found that a strong first-order QCD phase transition induces inhomogeneous 
nucleosynthesis. It is extremely important to understand the initial conditions for BBN, because many of 
our ideas about the early Universe rely on the validity of the standard (homogeneous) BBN scenario. This 
is in good agreement with observations 1351 . In inhomogeneous nucleosynthesis 1311 . large isothermal 
fluctuations of the baryon number (the remnants of the quark droplets at the end of the QCD transition) 
could lead to different yields of light elements. As a minimal requirement for an inhomogeneous scenario 
of nucleosynthesis, the mean bubble nucleation distance has to be larger than the proton diffusion length, 
which corresponds to ^ 3 m 1 32 1 at the QCD transition. This is two orders of magnitude above recent 
estimates of the typical nucleation distance [28 ). 

On the other hand the observed cosmic abundances of light elements do not favour inhomogeneous 
nucleosynthesis, except a small region in parameter space corresponding to an inhomogeneity scale of 
~ 40 m |32|. 

However, interesting inhomogeneity scales for BBN might follow if the bubble nucleation is not ho- 
mogeneous. Heterogeneous nucleation in the presence of impurities or nucleation in the inhomogeneous 
Universe could increase the baryon inhomogeneity scale to interesting values. 

Although values for 77 dramatically different from those in the standard BBN are excluded both from 
measurements of the light element abundances and from the CMB, it might be possible to alter the primor- 
dial abundance of heavy elements (A > 7) very much in various inhomogeneous scenarios li41J . More on 
inhomogeneous BBN will be discussed in section^ 



2.3.3 Cold dark matter clumps 

Scales A that are of the order of the Hubble radius i?H are not sensitive to details of the bubbles. It was 
reported in Refs. I50ll51ll52l that the evolution of cosmological density perturbations is strongly affected 
by a first-order QCD transition for subhorizon scales, A < i?H- Cosmological perturbations on all scales 
are predicted by inflation |87 88 89 90 91 1, observed in the temperature fluctuations of the cosmic 
microwave background radiation, first by the COBE satellite 1921 1931 . 

In the radiation-dominated Universe subhorizon density perturbations perform acoustic oscillations. 

The restoring force is provided by pressure gradients. These, and therefore the isentropic speed of sound 

1/2 

Cs — {dp/de)g (on scales much larger than the bubble separation scale) drop to zero at a first-order QCD 
transition 1511 . because both phases coexist at the pressure p^, only (a is the scale factor of the Universe): 

de(a)/da 

It stays zero during the entire transition and suddenly rises back to the radiation value Cs — 1/ \/3 after the 
transition. A significant decrease in the effective speed of sound Cs during the cosmological QCD transition 
was also pointed out by Jedamzik f491. 

As the speed of sound drops to zero, the restoring force for acoustic oscillations vanishes and density 
perturbations for subhorizon modes fall freely. The fluid velocity stays constant during this free fall. Per- 
turbations of shorter wavelengths have higher velocities at the beginning of the transition, and thus grow 
proportional to the wave number k during the phase transition. The primordial Harrison-Zel'dovich spec- 
trum 1 94 1 of density perturbations is amplified on subhorizon scales. The spectrum of density perturbations 
on superhorizon scales, A > i?H, is unaffected. At T ^ 1 MeV the neutrinos decouple from the radia- 
tion fluid. During this decoupling the large peaks in the radiation spectrum are wiped out by collisional 
damping |95!|. 

Today a major component of the Universe is dark matter, most likely CDM. If CDM is kinetically 
decoupled from the radiation fluid at the QCD transition, the density perturbations in CDM do not suf- 
fer from the neutrino damping. This is the case for primordial black holes or axions, but not for su- 
persymmetric dark matter. At the time of the QCD transition the energy density of CDM is small, i.e. 
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ecdm(7*) ^ 10^^ei-ad(7*). CDM falls into the potential wells provided by the dominant radiation fluid. 
Thus, the CDM spectrum is amplified on subhorizon scales. The peaks in the CDM spectrum go non- 
linear shortly after radiation-matter equality. This leads to the formation of CDM clumps with mass 
< lO^^°Af0. Especially the clumping of axions has important implications for axion searches 1961 . If the 
QCD transition is strong enough, these clumps could be detected by gravitational femtolensing 1971 . 

2.3.4 Other effects 

Generation of magnetic fields 

Hogan 1 42 1 argued that magnetic fields might be generated in violent processes at the bubble scale. 
Later on field lines should reconnect to yield large scale magnetic fields which could play an important 
role in structure formation after the recombination of matter. Cheng and Olinto |43| suggested that such 
magnetic fields might be generated by currents on the bubble surfaces. 

If the shrinking QGP droplets are baryon-enriched, there is a positive net charge on the inner side of the 
surface of the bubble walls, which is compensated by a negative net charge of electrons outside the bubble. 
This is due to the different Debye screening lengths of electrons and baryons. Thus surface cuiTents on the 
bubble surfaces are possible, which could give rise to magnetic fields. The typical scale at their creation 
would be the bubble scale. 

Subhorizon (A < Rn) fields are damped during neutrino and photon decoupling |98'i. In the most 
optimistic scenario, magnetic fields of the order of 10^^" gauB at the 10 Mpc scale would be possible 
today |99|. 

It might be that the QCD transition leads to the generation of pion-strings, predicted by an effective 
description of hadronic matter within the linear sigma model. These defects are unstable and decay even- 
tually, but could seed magnetic fields flOOl . The source of the magnetic field would come from the Adler- 
Bell-Jackiw anomaly, which couples the tt*^ to photons. 

A quite different mechanism has been proposed in IIOII . It relies on several untested assumptions, 
especially on the existance of axion domain walls and an inverse cascade mechanism to generate magnetic 
fields on scales much larger than 1 pc today. 

Generation of gravitational waves 

A measure for the energy density in gravitational waves is the fractional energy density per logarithmic 
frequency interval rigw(/) [see eq. ( I4.5l l1: / is the frequency of the gravitational wave under consideration. 

The generation of gravitational waves from violent bubble collisions during a strongly first-order QCD 
transition has been suggested by Witten |20|. The production of gravitational waves is possible if deto- 
nation bubbles collided, i.e. the bubble walls move faster than the speed of sound. A more quantitative 
analysis for the collision of vacuum bubbles (all latent heat goes into the kinetic energy of the bubble 
walls) has been performed by Kosovsky et al. 1 102 1. However, neither in the bag model nor from lattice 
QCD it is likely that a bubble nucleation scenario with detonation waves (almost vacuum bubbles) takes 
place. In the opposite, the most likely scenario is the occurence of weakly deflagrating bubbles (27|. The 
generation of gravitational waves under these circumstances has been investigated by Hogan 1 103| and 
Kamionkowski et al. 1 104|. An estimate in 1 103 1 yields figw ~ 10^^ {dnucH)^ c% for / ~ at the time 
of the QCD transition, which is / 10^^ Hz today. These gravitational waves are induced by the large 
inhomogeneities in energy density due to the coexistence of the QGP and the HG during a period ^ . 
The most optimistic scenarios allow ilg^ ~ 10~^^; however, scenarios based on recent lattice QCD results 
give figw ^ lO^'^'^, which is completely out of reach of any technique known today for the detection of 
gravitational waves. The calculations in 11041 result in a different dependence on dnucH; however, 17gw is 
as small as from 1 103|. 

Formation of black holes 

Crawford and Schramm 1471 suggested that long-range forces during the QCD transition lead to the 
formation of planetary mass black holes. However, there is no evidence that there are such large correlation 
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lengths at the QCD transition (km instead of fm!). Hall and Hsu 1481 argued that collapsing bubbles might 
collapse to black holes. The mass of these objects would be at most /oVbubbic, which, if the typical bubble 
scale is 1 cm, is < IO^^^AIq. Since the supercooling in the cosmological QCD transition is tiny, such a 
violent collapse of the walls of the shrinking quark droplets is impossible. Moreover, black holes of these 
masses could not have survived until today, they would have evaporated long ago fl05|. 

The vanishing of the speed of sound during the coexistence phase also leads to interesting gravitational 
effects, as pointed out by Jedamzik 1491 1 1061 . It was shown by several groups that the drop of the speed of 
sound is not sufficient to give rise to black hole formation 1521 1 1071 1 1081 . 



3 The radiation fluid at the QCD scale 

The expansion of the Universe is very slow with respect to the strong, electromagnetic, and weak interac- 
tions around (see figure [l~2i . To be more explicit, the rate of the weak interactions is ^ GpT^ « 
lO"^"* GeV, the rate of the electromagnetic interactions is Fcm ^ cP'Ti, w 10^^ GeV, and the rate of the 
strong interactions is Fg ^ a1{Ti,)Ti, « 10"^ GeV. These rates have to be compared with the Hubble 
rate H ~ /mp\ « 10~^^ GeV. Thus, leptons, photons, and the QGP/HG are in thermal and chemical 
equilibrium at cosmological time scales. All components have the same temperature locally, i.e. smeared 
over scales A ~ 10~'^i?H. At scales A > 10^^i?H, strongly, weakly, and electromagnetically interacting 
matter makes up a single perfect (i.e. dissipationless) radiation fluid. 

There are no conserved quantum numbers for the radiation fluid, apart from the lepton numbers and the 
baryon number However, the corresponding chemical potentials are negligibly small at the QCD epoch. 
For the baryon chemical potential this is shown below ("section 13741 1. All lepton chemical potentials are 
assumed to vanish exactly. In this situation all thermodynamic quantities follow from the free energy. The 
free energy density is f{T) — —p{T) from homogeneity. The entropy density is given by the Maxwell 
relation for the free energy: 

and the energy density e{T) follows from the second law of thermodynamics for reversible processes: 

e = T^-p. (3.2) 



3 . 1 Equation of state 

The behaviour of e(T) and p{T) near the QCD transition must be given by non-perturbative methods, by 
lattice QCD. In figure ITTI lattice QCD data are shown for e{T) (denoted by p in the figure) and p{T) 
divided by e of the corresponding SB gas. The lattice results for two systems are plotted: quenched QCD 
(no quarks) 1611 . and two-flavour QCD |70|. For quenched QCD the lattice continuum limit is shown. For 
two-flavour QCD the data with six time steps {Nf = 6, a « 0.2 fm) and a quark mass am^ ~ 0.0125 
is shown. This corresponds to a physical mass niq ~ 14 MeV, a bit heavier than the physical masses 
of the up and down quarks. On the horizontal axis we plot [T /T^,). For T/T^, = 4, energy density and 
pressure for quenched QCD are still 10% resp. 15% below the SB gas value. This is in excellent agreement 
with analytic calculations at finite temperature [110 111 |. It is remarkable that e/esB and p/esB versus 
T/Ti, is quite similar for quenched QCD and two-flavour QCD. Moreover, the temperature dependence of 
the rescaled pressure for four-flavour QCD |751 is quite similar to quenched QCD. For more lattice QCD 
results see |63|. At temperatures below quarks and gluons are confined to hadrons, mostly pions. At 
present the hot pion phase is not seen in the two-flavour lattice QCD, since the pion comes out too heavy 
(0.3 < m-^/mp < 0.7 from 1701 . while the physical ratio is 0.18). 
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Fig. 3.1 The energy density and the pressure of hot QCD, both relative to the energy density of a Stefan-Boltzmann 
(ideal) gas of quarks and gluons, are plotted as functions of T/T* (from |^). The lattice QCD data for Nf = b6ll 
and Nf — 2 1701 are compared with the predictions of the bag model. 



In the bag model (see section l2'.2.1> the pressure was given in eqs. ( I2.10> and (I2.11> and thus from (13. 1> 
and (13. 2> the energy density and entropy density read: 

eqcp = Eqcp + ^ (3-3) 

SQGP = 4gP' (3-4) 

with e^^ — {tt'^ /30)gT'^ and s^^ — {2Tr^ /45)gT^ for massless particles. Note that s/T^ is a simple step 
function in the bag model. 

3.2 Adiabatic expansion 

Entropy is conserved, apart from the very short stage of reheating (~ lO^^tn) after the first bubbles have 
been nucleated. This allows us to calculate T = T{a) from d{s{T)a'^) = 0, i.e. 

dT 3s 



din a ds/dr ' 



(3.5) 



except for T = in the case of a first-order phase transition. In the bag model, T cx 1/a for T ^T^. 

The expansion while the QGP and HG coexist in a first-order QCD transition is determined by entropy 
conservation, 

a+ = ' (3.6) 

where the index — (+) denotes the value of a quantity at the beginning (end) of the coexistence epoch. In 
the bag model the Universe expands by a factor of a+/a_ w 1.4 until all QGP has been converted into 
the HG, whereas for a lattice QCD fit |51 52| the Universe expands by a factor of a+/a_ « 1.1. The 
growth of the scale factor is related to a lapse in cosmic time by d In a = dt. In terms of the Hubble time 
the transition lasts 0.3tH, resp. O.I^h for the bag model resp. the lattice QCD fit of f5T"52^. 

During a first-order QCD transition, i.e. T ~ T*, the pressure p(r*) = is constant. For any first- 
order QCD phase transition the energy density p{a) is obtained from the first law of thermodynamics 
de = — 3(e + pj,)da/a. 
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Fig. 3.2 The behaviour of the speed of sound Cs — {dp/de}l^^ during the QCD transition as a function of a (from 
l52j ). During a first-order transition (lattice fit and bag model) the speed of sound vanishes. 



3.3 Speed of sound 

The speed of sound relates pressure gradients to density gradients, i.e. Vp = Ve. This relation is 
essential for the evolution of density fluctuations. During the short period of supercooling the relation 
between temperature and time strongly depends on c^. This is important for the correct estimate of the 
mean bubble nucleation distance. 

When analysing cosmological perturbations during the QCD transition, for wavelengths A > 10^'*i?H, 
neutrinos are tightly coupled: T^/k ^ 1. For these wavelengths the radiation fluid behaves as a perfect 
(i.e. dissipationless) fluid, entropy in a comoving volume is conserved, and the process is thus reversible. 
On the other hand, below the neutrino diffusion scale, A < acoustic oscillations are damped 

away before the QCD transition (see section l63i . 

The isentropic speed of sound (for wavelengths A much larger than the bubble separation), given by 

2 f 9p\ dp/da 

must be zero during a first-order phase transition for a fluid with negligible chemical potential (i.e. no 
relevant conserved quantum number), because p{a) = = constant. 

In the bag model, = 1/3 before and after the transition and vanishes during the transition^. For a 
crossover the speed of sound decreases below l\/3 at T^, but does not drop to zero (see 11 121 for a simple 
analytic model of a crossover transition). 

The isentropic condition applies during the part of the phase transition after the initial supercooUng, 
bubble nucleation, and sudden reheating to T^,. During the adiabatic part of the transition, which takes 
about 99% of the transition time, the fluid is extremely close to thermal equilibrium, because the time to 
reach equilibrium is very much shorter than a Hubble time, i.e. the fluid makes a reversible transformation. 
This can be seen as follows: across the bubble walls, local pressure equilibrium is established immediately, 
Pqgp = Pug locally. The local temperature equilibrium, Tqgp = Thg. is established by neutrinos, which 
have a mean free path of 10^^i?H, enormously larger than the bubble wall thickness, and a collision time 
much shorter than the Hubble time. This local pressure and temperature equilibrium can only be satisfied 



^ A small drop of the speed of sound was found by Dixit and Lodenquai 11131 in a bag model taking interactions and masses into 
account. However, they missed the fact that the speed of sound drops to zero in a cosmological first-order QCD transition. 
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if p = Pi, and T = T^, at the bubble walls. Over distance scales of the order of the bubble separation 
(~ 1 cm), pressure (and therefore also temperature) is equalized with the velocity of sound, and thereby 
the released latent heat is distributed. This pressure equalization is very fast with respect to the Hubble 
expansion velocity wh ~ lO^^c at the 1 cm scale. 

From eqs. (13. 2> and (13. 1> the speed of sound may be evaluated as 

= (dlns/dlnT)-i . (3.8) 

The relation between temperature and time during the adiabatic expansion depends on the speed of sound. 
From eqs. (13. 5> and (13. 8> we find 

dlnT = -SclHdt . (3.9) 

This relation holds during the short period of supercooling before the first bubbles have been nucleated. 
An estimate of in the QGP at T^, has been obtained from lattice QCD. 

A strong decrease in the speed of sound, already above T^, has been observed in lattice QCD for Nf ~ 
161) and for Nf = 2 II70II66I . In both cases, Cg(Tt) w 0.1 when approaching the critical temperature from 
above (in the QGP). Note that this implies that the speed of sound relevant to bubble nucleation is of order 
O.Sinsteadof w 0.6. 



3.4 Baryons 

The baryons are tightly coupled to the radiation fluid at the QCD scale. Their energy density is negligible 
with respect to that of the relativistic particles (photons, leptons, quarks/pions), thus they are dragged with 
the radiation fluid. Below 1 argue that the baryon chemical potential is negligible at the QCD scale. 

At high temperatures (T > T^) the baryon number density may be defined as ne = ^ ~ ''M)' 

where nq (riq) is the number density of a specific quark (antiquark) flavour, and the sum is taken over all 
quark flavours. At T < 1 GeV only the u, d, and s quarks contribute significantly. At low temperatures 
(T < Ti,) the baryon number density is defined as n-g, = X](^*b — n^^)', now the summation is taken over 
all baryon species. Practically the nucleons contribute to the baryon number of the Universe only. 

Below the electroweak transition (Tew = 100-200 GeV) the baryon number i? in a comoving volume 
is conserved. On the other hand, the entropy S is conserved. As a consequence the ratio of baryon number 
density and entropy density s is constant. From the abundances of primordial ^He and D we know the 
ratio TiB /n-y — (5.6 ± 0.5) x 10~^" |38|. Taking into account the three massless neutrinos along with the 
photons that contribute to the entropy density, we find 

— = (7.5±0.7) X 10"" . (3.10) 
s 

Owing to the smallness of this ratio, the number of quarks equals the number of antiquarks in the very 
early Universe. 

Let me now turn to the baryon chemical potential. At high temperatures the quark chemical potentials 
/iq are equal, because weak interactions keep them in chemical equilibrium (e.g. u + e<-^dors + ly^), and 
the chemical potentials for the leptons are assumed to vanish (see 1 109 1 for a discussion of lepton chemical 
potentials). Thus, the chemical potential for a baryon is defined by /ie = 3/iq. For an antibaryon the 
chemical potential is —pB- The baryon number density of an SB Fermi gas of three quark flavours reads 
riB ~ T^/xb/S at high T. From eq. (I3.10> one finds that 

- 10"^ atr>r,^. (3.11) 

At low temperatures (T < T*), ij,b — fJ-p — l^-n, neglecting the mass difference between the proton and 
the neutron. The ratio of baryon number density and entropy now reads 




(3.12) 
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Since this ratio is constant, the behaviour of ^-b/T is given by eq. 
Hb/T « 10^^. All antibaryons are annihilated when the ratio 



e.g. at iTip/T « 20 one finds 



nh - »b 



tanh ( ^) 



(3.13) 



goes to unity. This happens when ijl-b/T ~ 1, which corresponds to Tann ~ fTip/25 « 40 MeV. Below this 
temperature the baryon chemical potential is ijlb{T ^ Tann) ~ 'rrip- To add one proton to the Universe 
one proton rest mass should be invested. A detailed investigation of the baryon chemical potential in the 
early Universe was recently given by 111 141 . 

It was argued above that, from sphaleron processes before the electroweak transition, a possible lepton 
asymmetry should be of the same order of magnitude as the baryon asymmetry. Thus, at temperatures 
T ^ me we find that lepton chemical potentials are negligible, since ^i^/T ^ ni^/T^ ^ ub/T^ ^ t]. 
This justifies the previous assumption that any chemical potential can be neglected during the first second. 



4 Evolution of gravitational waves 

In principle, primordial gravitational waves (e.g. from cosmological inflation) present a very clean probe 
of the dynamics of the early Universe, since they know only about the Hubble expansion. As was shown 
in 1531 a step is imprinted in the spectrum of primordial gravitational waves by the cosmological QCD 
transition. This step does not allow us to tell the difference between a first-order transition and a crossover, 
but its position would allow an estimate of the temperature and its height would allow a measurement of 
the change in the effective number of relativistic degrees of freedom. 

Primordial gravitational waves are predicted to be generated during inflation fl 15' '1 161 and could be 
detected by observing the so-called B-mode (parity odd patterns) polarisation of the CMB. Inflation pre- 
dicts an almost scale-invariant energy density per logarithmic frequency interval for the most interesting 
frequencies (- IQ-^ Hz for pulsar timing, - lO^^ Hz for LISA, and - 100 Hz for LIGO and VIRGO) of 
the gravitational waves. 

In the cosmological context, the line element of gravitational waves is given by 

ds^ = -dt^ + a2(5,y + h,j)dx'dx^ , (4.1) 

where hij is a transverse, traceless tensor The linearised Einstein equation admits wavelike solutions for 
hij^. The spatial average 

{h,j{x)h'^x + r)) = J ja{kr)k^\hk\^dlnk (4.2) 

defines the power spectrum |/ifcp. We denote by h the rms amplitude of a gravitational wave per logarith- 
mic frequency interval: h = k'^/^lhk]. The linearized equation of motion for h{t) reads 

1.2 

h + '3Hh +—h^0, (4.3) 

where the differentiation is taken with respect to cosmic time t. The amplitude of gravitational waves is 
constant on superhorizon scales and decays as 1/a after horizon crossing, h ~ Ck sin(fcr7 + dk)/a, where 
77 = J dt/ais conformal time; Ck and Sk are determined by matching the subhorizon to the superhorizon 
solution. 



For a more general definition and discussion of gravitational waves see Ref. ITTTl . 
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Fig. 4.1 The modification of the energy density, per logarithmic frequency interval, for primordial gravitational waves 
from the QCD transition (from 1531). 



For subhorizon modes, fcph ^ H, the energy density of gravitational waves can be defined. The space- 
time average of the energy-momentum tensor over several wavelengths gives egw — ~{l/32TTG){hijh^^). 
The energy density per logarithmic interval in k is related to the rms amplitude h: 



, dcgw 



327rG P''2 



(4.4) 



The factor 1/2 comes from the time average over several oscillations. The energy fraction in gravitational 
waves, per logarithmic interval in k, is defined by 

degw 1 



^gw(^) — ^" 



dfc 



(4.5) 



where = SH^/SnG. 

Figure ini shows the transfer function f^gw(/)/^^gw(/ ^ /*) from a numerical integration of eq. ( 14. 3> 
through the cosmological QCD transition. The typical frequency scale is 



A « 1.36 



9 



rlO^^Hz, 



(4.6) 



.17.25/ 150 MeV 

which corresponds to the mode that crosses the Hubble horizon at the end of the bag model QCD transition. 
Scales that cross into the horizon after the transition (l.h.s. of the figure) are unaffected, whereas modes 
that cross the horizon before the transition are damped by an additional factor w 0.7. The modification of 
the differential spectrum has been calculated for a first-order (bag model) and a crossover QCD transition. 
In both cases the step extends over one decade in frequency. The detailed form of the step is almost 
independent from the order of the transition. 

The size of the step can be calculated analytically |53 1. Comparing the differential energy spectrum for 
modes that cross into the horizon before and after the transition gives the ratio 



^^gw(/<A) V.9b 



0.696 , 



(4.7) 
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for the QCD transition, which coincides with the numerical integration in fig. l4.1l This resuh is in agree- 
ment with the entropy conservation of subhorizon gravitational waves I103II1 181 171. However, for super- 
horizon modes the entropy is not defined. 

Similar steps in the differential spectrum have been studied for gravitational waves generated by cosmic 
strings il.19 1. These gravitational waves are generated on subhorizon scales when the cosmic strings decay. 
Their frequency is larger than a^^H, where a is the ratio between the typical size of a string loop and the 
Hubble radius at formation of the loop; a is smaller than 0.1 and might be as small as 10^'' 1 120|. In 
this situation, steps in the differential spectrum follow from the conservation of entropy for decoupled 
species (gravitons) during cosmological phase transitions L121J . This interpretation applies to modes that 
have been inside the horizon long before the transition or that have been generated on subhorizon scales. 
However, for superhorizon modes, entropy and energy of a gravitational wave are not defined. We therefore 
cannot rely on a conservation-of-entropy argument when dealing with gravitational waves from inflation. 

In fig. O we indicated the frequency range 1 yr ^) in which limits on fig^ have been reported 
from pulsar timing residuals 1 122 1. The frequencies where the step of the QCD transition would be visible 
is of the order of 0.3 month^^. For pulsar timing, the power spectrum of gravitational waves is more 
relevant than the energy spectrum. The power spectrum is cx ilgw(/)/^'^. Our results show that the power 
spectrum might deviate from the /^^ behaviour over a whole decade in frequency. Unfortunately, with 
today's technology we will not be able to detect primordial gravitational waves at frequencies around 10^^ 
Hz, because their amplitude is expected to be to small. 



5 A first-order QCD transition 

In a first-order QCD transition the quark-gluon plasma supercools before the first bubbles of hadron gas 
are formed. In a homogeneous Universe without 'dirt' the bubbles nucleate owing to thermal fluctuations 
(homogeneous nucleation). 

If cosmic 'dirt' in the form of defects (such as strings) or black holes is present at the QCD epoch, this 
'dirt' may trigger the formation of the first hadronic bubbles (like the nucleation of vapour bubbles in a pot 
of boiling water). 

There is a broad range in parameter space, where the magnitude of primordial temperature fluctuations 
is of the same order or larger than the typical supercooling in the homogeneous nucleation scenario. In this 
case, the transition proceeds with inhomogeneous bubble nucleation. The mean nucleation distance results 
from the scale and amplitude of the temperature fluctuations. 



5 . 1 Homogeneous nucleation 

The probability to nucleate a bubble by a thermal fluctuation is proportional to exp(AS'), where AS is the 
change in entropy by creation of a hadronic bubble. The second law relates AS* to the minimal work done 
in this process, which is the change in the free energy because the volume and temperature are fixed. The 
change in free energy of the system by creating a spherical bubble with radius R is 

AF = —{PQGP - Phg)-R' + 47rai?2 , (5.1) 

where a is the surface tension. Bubbles can grow if they are created with radii greater than the critical 
bubble radius Re- Smaller bubbles disappear again, because the free energy gained from the bulk of the 
bubble is more than compensated by the surface energy in the bubble wall; Rc is determined from the 
maximum value of AF{R) and reads 

Rc{T) = r-T^ r-: ■ (5-2) 
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At Ti, the critical bubble size diverges, and no bubble can be formed. Finally, the probability to form a 
hadronic bubble with critical radius per unit volume and unit time is given by 

/(r)=/o(r)exp(^-^^ , (5.3) 

with AFc = 167r(T'^/[3(j3HG^J5QGp)^]. For dimensional reasons the prefactor /q ^ CTf , with C = C(l). 
A more detailed calculation of /q within the bag model has been provided in |123|. It was shown in 
Ref. 1 28 1 that the temperature dependence of the prefactor Jg can be neglected for the calculation of the 
supercooling temperature Tgc in the cosmological QCD transition. As will be shown below, the numerical 
prefactor C is iiTelevant in the cosmological QCD transition. 

For small supercooling A = 1 — T/T* ^ 1 we may evaluate (phg ~ PciGp){T) by using the second 
law of thermodynamics, i.e. phg ^ Pqgp ~ and thus 

/(A)«/o(TOexp(-A/A2), (5.4) 

with A = 167ro-3/(3Pr^) and Iq{T^,) k, Tf. Note that this result does not depend on the details of the 
QCD equation of state. For the values of I = 1.4T^ and a = O.OIST^ from quenched lattice QCD |25| 
^ « 3 X IQ-^. In the bag model A « 5 x 10~2(cr/7;3)3. 

The amount of supercooling that is necessary to complete the transition. Age, can be estimated from the 
schematic case of one single bubble nucleated per Hubble volume per Hubble time, which is 



0(As 



A 



A\u{T./H,) 



1/2 



4 X 10~^ (5.5) 



for the values of I and a from quenched lattice QCD |25|. For the bag model I assume a < O.IT^, which 
implies that Age < 6 x 10^''. It has been shown in |27| that for such a tiny supercooling the formation of 
detonation bubbles is forbidden. This justifies the approximation of small supercooling made above. 
The time lapse during the supercooling period follows from the conservation of entropy and reads 

Atgc/tH = Asc/(3c2) = 0(10-3) . (5.6) 

Here we used the relation = d In s /d In T for the speed of sound in the supercooled phase. For realistic 
models < Cs(A) < l/VS. In the bag model Cs( A) = l/\/3. 

The critical size of the bubbles created at the supercooling temperature is 

i?c(A,e)«|^«25fm. (5.7) 

The critical radius is large with respect to the QCD scale, and this justifies the thin wall approximation, 
which was made implicitly above. 

After the first bubbles have been nucleated, they grow most probably by weak deflagration I19ll23ll24l 
|271. The deflagration front (the bubble wall) moves with the velocity VdcH <C l/\/3 [1241. The energy 
that is released from the bubbles is distributed into the suiTounding QGP by a supersonic shock wave 
and by neutrino radiation. This reheats the QGP to and prohibits further bubble formation. Since the 
amplitude of the shock is very small f231, on scales smaller than the neutrino mean free path, heat transport 
by neutrinos is the most efficient. Neutrinos have a mean free path of 10-^i?H at T*. When they do most 
of the heat transport, heat goes with t^heat = C(c). For larger scales, heat transport is much slower. Figure 
l5.1l shows a sketch of the homogeneous bubble nucleation scenario. 

Let us now calculate the mean bubble separation, c?nuo and the final supercooling. Age, for a scenario 
with weak deflagration. Bubbles present at a given time have typically been nucleated during the preceding 
time interval 

AW = //(d//dt) . (5.8) 
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Fig. 5.1 Sketch of a first-order QCD transition via homogeneous bubble nucleation: above the critical temperature 
the Universe is filled with a quark-gluon plasma (Q). After a small amount of supercooling the first hadronic bubbles 
(H) nucleate at t\, with mean separation dnuc- At > t\ these bubbles have grown and have released enough latent 
heat to quench the formation of new bubbles. The supercooling, bubble nucleation, and quenching takes just 1% of the 
full transition time. In the remaining 99% of the transition time the bubbles grow following the adiabatic expansion 
of the Universe. At the transition is almost finished. The shrinking QGP drops are separated by the typical distance 



Using the relation between time and supercooling, dA/di — SCg/tn, we find 

AW/te = A3j(6Ac2) = 0(10-^) (5.9) 

and 

A„„c = ^^Asc = O(10-2)A,^ . (5.10) 

During the time interval A<nuc each bubble releases latent heat, which is distributed over a typical distance 
« A2whcatinuc- This distance has a weak dependence on the precise value of Age, but the bubble nucleation 
rate increases strongly with A until one bubble per volume ^ (Ainuc^'hcat)'^ is nucleated. Therefore the 
mean bubble separation is 

dnuc ~ 2t;heatAW « ^^^H = 0{\^-^ R^i) = O(lcm), (5.11) 

where I used Wheat = C(O.l), 3c^ = 0(0.1), which gives a typical value for the nucleation distance. The 
suppression of bubble nucleation due to already existing bubbles is neglected. 

The estimate ( 15.1 1> of the mean bubble separation applies if the released latent heat by means of sound 
waves and by neutrino free streaming is sufficient to reheat the QGP to T*, i.e. to quench the nucleation 
of new bubbles. On the other hand the typical bubble separation could be given by the rate of release 
of latent heat, i.e. by the bubble wall velocity Wdofl- Since the period of supercooling lasts about 1% of 
the time needed for completing the entire first-order phase transition, 1% of the QGP must be converted 
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to HG in the process of sudden reheating to T^; the bubble radius at quenching must therefore reach 0.2 
of the bubble separation, i?bubbic ~ 0.2(i„uc- With i?bubbic ~ WdcflAinuc, and using the above relation 
dnuc ~ 2whcat Aiiiuc, we require Udofi > 0.4uhcat for consistency. If Vdcfi is smaller than this, the limiting 
factor for quenching is the rate of release of latent heat by bubble growth, and the bubble separation is 

dnuc ~ 27;dcfiAi„uc « ^^^^i?H , (5.12) 

i.e. the bubble separation will be smaller than the estimate in eq. ( I5.11> . 

We are now in a position to improve the estimate of Age : one bubble nucleates in the volume (whcat Atnuc ) 
during Atnuc- This can be written as 

1 « (i;hcatAt„uc)^AW/(tsc) , (5.13) 
which in terms of the supercooling parameter Age is given by: 



3 



lO'^'A'J exp ( i . (5.14) 



X 10 



-5 



Also the pre-exponential factor is smaller by a factor of 10^'^ than the naive estimate i5.5i . the amount of 
supercooling is just 20% larger than in i5.5\ . i.e. Age = 5 x 10^^. This also demonstrates that numerical 
prefactors in ( I5.13> are irrelevant in the calculation of Age. 

To summarize, the scales on which non-equilibrium phenomena occur are given by the mean bubble sep- 
aration, which is about 10^^i?H- The entropy production is tiny, i.e. AS/S ~ 10~®, since the supercooling 
is small ~ 10~^. After supercooling, which lasts lO^'^tn, the Universe reheats in Atnuc ~ lO^^in- After 
reheating, the thermodynamic variables follow their equilibrium values and bubbles grow only because of 
the expansion of the Universe. 



5.2 Heterogeneous nucleation 

In the first-order phase transitions that we know from our everyday experience, for example the conden- 
sation of water drops in clouds, the drops, i.e. the bubbles are nucleated at impurities ('dirt'). This could 
happen in the early Universe as well. Candidates for cosmic 'dirt' are primordial black holes, monopoles, 
strings, and other kinds of defects. Of course, the existence of any of these objects has not been verified 
so far Nevertheless, let me discuss in what manner cosmic 'dirt' would change the nucleation of bubbles. 
The following considerations are based on the work of Christiansen and Madsen 1281 . 

Let n be the number density of the impurities. Further assume that at time ti with t-^ < ti < tgc bubbles 
nucleate at the locations of the impurities. (It is easy to see that is restricted to the mentioned interval: 
before ti, a bubble cannot grow because T > T*, after tgc homogeneous nucleation already happened.) 

There are two limiting cases: if n ^ d^^^ the mean nucleation distance is dnuc = n^^^^ 

dnuc, horn- If ?T- < d~^^^ the impurities are so rare that d^uc = rfnuc, hom- 

The most interesting situation occurs when the typical distance between the impurities is something big- 
ger than the mean homogeneous nucleation distance. But it should be small enough for the bubbles nucle- 
ated at the impurities to reheat the Universe just before homogeneous nucleation starts. For a quantitative 
estimate, let me determine the amount of supercooling for heterogeneous nucleation. As in homogeneous 
nucleation [see eq. ( I5.13H I use the condition that one bubble is nucleated per reheated volume, i.e. the 
sum of the probabilities to form a bubble from an impurity and from a thermal fluctuation: 

1 ~ n(whcatAt)3 + (t>hcatAtnuc)^At„uc/(tsc) ■ (5.15) 
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At — tsc — ti is the growth time of bubbles from impurities and Atnuc the growth time of bubbles from 
thermal fluctuations, as defined in ( 15. 8t . Note that here <sc may be smaller than in homogeneous nucleation. 
If At/Ai„uc is of order unity both mechanisms work independently and 

_ ^nuc, hom ■ 

(5.16) 



(n + AW/(isc))i/3 

The most interesting situation arises when At/At„uc ^ 1. This means that the probability for a bubble 
from a thermal fluctuation is reduced because the available volume for thermal fluctuations in the QGP has 
been reduced by a factor of (1 — n(whoat At)'^). This factor should be taken into account in (I5.15> which 
now reads 

1 ~ n(fhcatAt)3 + (whoatAW)^At„uc/(tsc)[l " n(vhoat At)^] . (5.17) 

The unique solution to this equation is n(whcat At)'^ « 1 and thus 

^ n-^l^ . (5.18) 



{n + AW/(tsc)[l - n(i;heatAt)3]}i/3 
The maximal mean nucleation distance from heterogeneous nucleation is found from At < Aigc. hot to be 



max(dnuc) ~ -yhcatmax(Atsc, hot) ~ WhcatAtgc, hom ~ 10 ^i?H = 0(10 m) . (5.19) 

There, for a narrow range in parameter space, namely n(i;hcat At)^ sa 1, nucleation distances maybe larger 
than in homogeneous nucleation by a factor of 100, which means that n should be n ~ 10^/ i?^. 

As an example for 'dirt' 1 consider primordial black holes (PBHs): with the appropriate density, 
'T.pBHlr^) ^ 10^/i?H, PBHs are produced at the temperature Tpbh ^ 10^ GeV [see eq. dO l with a 
mass of mpBH ^ lO~^®A/0. PBHs of this mass are excluded observationally because they start to evapo- 
rate today and should be observed as 7-sources [105 J. Thus PBHs cannot be considered as seeds that yield 
the maximum nucleation distance. Nevertheless, larger PBH masses still may seed bubbles, but their mean 
separation is so large that thermal bubble nucleations cannot be suppressed. 

Another example of 'dirt' are cosmic strings. Recently, such a scenario has been analysed in [125J . A 
moving cosmic string would generate an overdense plane, on which the phase transition would be delayed. 

5.3 Inhomogeneous nucleation 

The local temperature T(t,x) of the radiation fluid fluctuates, because cosmological perturbations have 
been generated during cosmological inflation |90|. Let me denote the temperature fluctuation by A^ = 
5T /T. Inflation predicts a Gaussian distribution of perturbations: 

1 / 1 A^ \ 

P(AT)dAT = ^= cxp — ^ cIAt . (5.20) 

^ ' \/2^A™- 2(A™^)2; 

If one allows for a tilt in the power spectrum of density fluctuations, the COBE 1931 normalized rms 
temperature fluctuation reads If29l 

/ L \ ("-l)/2 

A?i'-«10-4(3c2)3/M^j , (5.21) 

where feg is the wave number of the mode that crosses the Hubble radius today. The case n = 1 gives 
the Harrison-Zel'dovich spectrum |94|. Recent WMAP results combined with other CMB data and 2dF 
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Fig. 5.2 Sketch of a first-order QCD transition in the inhomogeneous Universe (from |29J): at ti the coldest spots 
(dark grey) are cold enough to render the nucleation of hadronic bubbles (H) possible, while most of the Universe 
remains in the quark-gluon phase (Q). At t2 > ti the bubbles from the cold spots have merged and have grown to 
bubbles as large as the fluctuation scale. Only the hot spots (light grey) are still in the QGP phase. At ts the transition 
is almost finished. The last QGP drops are found in the hottest spots of the Universe. The mean separation of these hot 
spots can be much larger than the homogeneous bubble nucleation separation. 



galaxy redshift survey results give n — 1 = —0.03 ± 0.03 (ti — 1 = — 0.07_q q5) without (with) running of 
the spectral index (dn/d In fc ^ 0) |3|. For n = 1 we find A'rj.''^ (kQco) ~ 2 x lO"'^. 

From the above we conclude that A™** and Anuc [see eq. ( 15.1 OH may be of the same magnitude or 
that A™^ may be even larger. The picture of homogeneous bubble nucleation, where bubbles form from 
statistical fluctuations, is false for the most probable cosmological scenarios. 

We thus developed a new scenario for the cosmological QCD transition |29|. To do that we had to 
learn more about the primordial temperature fluctuations first. A small scale cut-off in the spectrum of 
primordial temperature fluctuations comes from collisional damping by neutrinos I95II52I (see also section 
16. 3> . Because the neutrinos interact only weakly, their mean free path is large with respect to the strong 
and electromagnetic interacting particles. The interaction rate of neutrinos is F ~ G|T^. This has to be 
compared with the frequency Cgkph- We find that at the QCD transition neutrinos travel freely on scales 
Ai/-mfp ~ 10^^i?H- Fluctuations on the diffusion scale of neutrinos are washed out by the time of the 
QCD transition (see section l63t : 

1 _4 

Ajy-diff = g V ^i^-mfp'^^H ~ 10 -Rh • (5.22) 
Thus the old picture of homogeneous bubble nucleation still applies within the small homogeneous patches 

of Asmootti — 10~''i?H- 

The compression time scale for a homogeneous patch is 6t — Agmootii/cs ^ lO^'^in- If the compression 
time scale is larger than Atnuc the temperature fluctuations are frozen with respect to the time scale of 
nucleations. 

A sketch of inhomogeneous bubble nucleation is shown in fig. 15.21 The basic idea is that temperature 
inhomogeneities determine the location of bubble nucleation. In cold regions, bubbles nucleate first. In 
general we have two possible situations: 
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1. If A„uc > A™^, the temperature inhomogeneities are negligible and the phase transition proceeds 
via homogeneous nucleation (see section lsHl . 

2. If Anuc < A™^, the nucleation rate is inhomogeneous and we have to consider the scenario sketched 
in fig.lO 

A first attempt to analyse inhomogeneous nucleation has been given in f29V According to f29], the nucle- 
ation distance dnuc exceeds the scale Asmooth, if 

Asmooth < 2-|^A™^i?H- (5.23) 

If A™^ > 5 X 10^^, it is quite likely that this condition is met. In that case we can conclude that the typical 
inhomogeneity scale in the baryon distribution is inherited from the scale of density inhomogeneities in the 
radiation fluid at the end of the QCD transition. The effect in terms of length scales is at least two orders 
of magnitude larger than the nucleation distance in homogeneous nucleation and is 0(1 m), which is of 
interest for inhomogeneous BBN. 



6 Density fluctuations of the radiation fluid 

6. 1 Amplification of fluctuations 

Since the speed of sound vanishes during a first-order QCD transition (see section lSJt . the restoring forces 
for compressional perturbations vanish and thus density inhomogeneities on scales below the Hubble scale 
are amplified liSTl 1521 [108 1. For small perturbations e{x,t) =: eo(i) + Se{x,t), the equations of motion 
can be linearized in the perturbations. Here we are interested in the density perturbations, the quantity of 
interest is the density contrast 5 = i5e(x, t)/€o{t). 

The transfer functions, i.e. the change in the primordial spectrum, for the radiation fluid and the cold 
dark matter (CDM) are calculated in I5T1I52I. 

For the bag model the transfer functions are shown in fig. 16. II Both transfer functions show huge peaks 
on small scales. The different scales /cph = 27r/A are represented by the invariant CDM mass contained in 
a sphere of radius A/2, 



Mcdm(A) = 




(6.1) 



assuming for simplicity that fldm ~ 1 today. The largest scales in fig. 16. II correspond to the horizon 
at Tf = T*/10. The CDM curve also shows the logarithmic growth of subhorizon scales of CDM in a 
radiation-dominated Universe. The CDM peaks lie on top of this logarithmic curve. 

The peak structure starts at a scale « lO^^A/0 in CDM mass. This scale corresponds to the horizon 
scale at the QCD transition. The radiation energy inside the horizon at T* is ^ IMq, but it is redshifted as 
Mra.d{a) ~ (acq/a)A/cdm. Scales which are above the horizon at the QCD transition are not affected. For 
scales below M^'^'^ « 9 x IO^^A/q the radiation peaks grow linearly in wave number This linear growth 
comes from the fact that the vanishing speed of sound during the QCD transition implies a vanishing 
restoring force in the acoustic oscillations on subhorizon scales. Therefore, the radiation fluid falls freely 
during the transition, with a constant velocity given at the beginning of the transition. The density contrast 
^lad gj-Qws linearly in time with a slope k. CDM is moving in an external potential provided by the 
dominant radiation fluid, and is pushed by the strong increase in the gravitational force during the transition. 
The highest peaks have k/ki ^ 10^, because on smaller scales the acoustic oscillations are damped away 
by neutrino diffusion already before the QCD transition (see section l63l . 

The processed spectrum for a crossover, fig. 16.21 shows a behaviour similar to that for the bag model 
on superhorizon and horizon scales. The peak structure starts at AIi, but on subhorizon scales there are no 
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Fig. 6.1 The modifications of the density contrast for kinetically decoupled CDM (e.g. axions or primordial black 
holes), ^cdm = 15cdm|(r*/10), and of the radiation fluid amplitude, A^ad = (^rad + ^'ipLd)^^^' due to the QCD 
transition in the bag model (from l(52J ). Both quantities are normalized to the pure Harrison-Zel'dovich radiation 
amplitude. On the horizontal axis the wave number k is represented by the CDM mass contained in a sphere of radius 
Tv/k. 
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Fig. 6.2 The same as fig. 16. II but for a QCD crossover (from 1521 '). 

peaks. The level of the subhorizon transfer function for the radiation fluid is reduced to 0.83. This comes 
from the damping of the acoustic oscillations during the time with 1/3. 

The time evolution of subhorizon modes, A <C Rh, can be solved analytically during the transition. For 
the dynamics of the radiation fluid (QCD, photons, leptons) CDM can be neglected, since ecdm/eiad = 
a/oeq ~ 10~®. The transition time is short with respect to the Hubble time at the transition, (i+ — 
t-) < tn = H^^. For subhorizon modes we can neglect gravity during the whole transition, as has been 
shown in |52|. The damping terms in the continuity equation and Euler equation are absent in the purely 
radiation dominated regime. During the transition the damping terms can be neglected in view of the huge 
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amplification for a first-order phase transition. The continuity and Euler equations read: 

S' -k^p^O (6.2) 
ip' + clk5 = , 

where %!} denotes the peculiar velocity up to a factor e/{e + p). The primes denote derivatives with respect 
to conformal time. Written as a second-order differential equation for 5, this is just an oscillator equation 

S" + = , (6.3) 

with the time-dependent frequency = kcg{ri). 

Let us discuss the origin of large amplifications for small scales in the bag model. Before and after 
the QCD transition, the radiation fluid makes standing acoustic oscillations in each mode fc, with speed of 
sound (?g — 1/3 and with amplitudes A-^a and Aout for the incoming and outgoing solution respectively, 
see fig. 16.31 The incoming solution for the density contrast 5 and the peculiar velocity \/3V^ reads 

5 = - Ain cos [lu {ri - 7]-) + ip-] (6.4) 
V3-0 = AinSin[Lu{ri - T]^) + ip^] . 

This solution is valid until the beginning of the transition at rj = r/^, and denotes the phase of the 
oscillation at rj^. During the transition the speed of sound is zero. There are no restoring forces from 
pressure gradients and the radiation fluid falls freely. Since the duration of the transition is short with 
respect to the Hubble time At = (t+ — i_) < H^^, gravity is negligible during this free fall. The fluid is 
thus moving inertially in the sense of Newton, the velocity stays constant, and the density contrast grows 
linearly in time: 

S = 6- + k{r] - f]^)ip_ , (6.5) 

where VStp- = Ain sin (if-) is the peculiar velocity at rj-. 

Since we have no jumps in pressure, the density contrast S and the fluid velocity t/; stay continuous 
throughout the whole transition, in particular at the matching points of the different regimes. Gravity 
remains negligible during the entire transition. At the end of the transition this solution has to be joined 
to the pure radiation-dominated regime for T < T*. Since the ampUtude of the density contrast grows 
linearly during the transition, the final ampUtude ^out — A^ is enhanced linearly in k, modulated by the 
incoming phase 

^^ysin2(^_), (6.6) 



Ain J \ki 

with fci = VS/ A77, Arj = rj^ ~ -q^. The envelope of the linearly growing peak structure for subhorizon 
scales starts at the scale fci, which corresponds to a CDM mass of Mi = 9 x lO^^A/©. 

In the case of a crossover, the amplification occurs for scales around the Hubble radius at the transition 
only. Subhorizon scales always stay in the WKB-regime, and therefore the spectrum is flat for these scales. 
However, the amplitude for subhorizon scales is damped during the phase transition. The same damping 
occurs in the case of a first-order phase transition. It has been neglected in the analytic discussion, since it 
is a small correction. The subhorizon amplitudes are reduced to 83% of their initial value >52il . 

6.2 Formation of black holes? 

It was suggested in the literature 1471 1491 [T06l that the QCD transition could lead to the formation of \Mq 
black holes, which could account for today's dark matter Jedamzik 1491 proposed to identify such black 
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Fig. 6.3 The time evolution of the density contrast 5rad and the peculiar velocity l-^rad of the radiation fluid for the 
mode k/ki — 7 (from L521 ). During the QCD transition in the bag model — marked by the 2 vertical lines — the 
velocity stays approximately constant and the density contrast grows linearly. The amplitude is normalized to 1 long 
before the transition. 

holes with the MACHOs observed by microlensing 1861 . He pointed out that the formation of black holes 
should be particularly efficient during the QCD epoch thanks to the significant decrease in the effective 
speed of sound. 

In order to form a black hole in a radiation-dominated Universe, the density contrast inside the Hubble 
radius should be in the range 1/3 < (5h < 1 11261 . For an observable amount of IMq black holes 
today, i.e. Jlg^ = 0(1), the fraction of energy density converted to black holes at the QCD transition 
must be 0(aQCD/aoq) ~ 10^^. For a Gaussian distribution this requires 5,.,„s « 0.06 (without including 
any enhancement from the QCD transition) 1 127 1. The QCD transition gives an enhancement factor (at the 
horizon scale) of 2 for the bag model and of 1 .5 for lattice QCD in the linear perturbation treatment, figs. 16. II 
and 16. 21 This indicates a corresponding reduction in the required pre-existing perturbation spectrum at the 
solar mass scale. Cardall and Fuller used a qualitative argument of CaiT and Hawking 1 128 1 and the bag 
model; they also obtained a reduction of a factor of 2 in the required pre-existing perturbation spectrum 
1 107 1 . These QCD factors of 1.5 or 2 are so modest that a pre-existing Harrison-Zel'dovich spectrum with 
COBE normalization is very far from giving a cosmologically relevant amount of black holes II108I . One 
would have to put in a fine-tuned tilt (n — 1) « 0.4 to get the desired amount of black holes. However, this 
tilted spectrum would overproduce primordial black holes on scales that are only a factor of 50 below the 
Hubble radius at the QCD transition. Therefore a break in the pre-existing spectrum below the QCD scale 
would be required, a second fine-tuning. 

We conclude that the QCD transition does not produce black holes, although it enhances the probability 
to form some. The pre-existing spectrum needs to be fine-tuned around the QCD scale if one would like to 
produce a detectable amount of black holes. Thus the major effect would not be due to the QCD transition, 
but due to a feature in the primordial spectrum. 

6.3 CoUisional damping at neutrino decoupling 

The acoustic oscillations in the radiation fluid get damped by neutrino diffusion at the time of neutrino 
decoupling. This damping is analogous to Silk damping at photon decoupling. The muon and tau neutrinos, 
which are coupled to the radiation fluid via neutral current interactions only, decouple at T^^^^^ ~ 2.2 MeV 
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from the Hubble scale i?H, which follows from Ref. 1 129|. The electron neutrinos interact by charged and 
neutral currents and decouple slightly later, T^"'^ ~ 1.4 MeV. By the time of neutrino decoupling at the 
Hubble scale all inhomogeneities in the radiation fluid on scales below sa IO^^AIq in CDM mass are 
wiped out by neutrino diffusion (cf. the QCD horizon scale is 1Q~^Mq in CDM mass), as shown below, 
eq. ( 16. lU . Therefore the QCD peaks cannot affect big bang nucleosynthesis (BBN). 

It is important to distinguish the total decoupling of neutrinos, i.e. a neutrino decoupling at the Hubble 
scale, when the neutrino scattering rate T is smaller than the Hubble rate, T^/ H < 1, from the neutrino 
decoupling with respect to a certain mode given by T^/ujpi^ < 1, when the typical neutrino scatters less 
than once during an acoustic oscillation time of one particular mode. The mode-dependent decoupling 
temperature T^^^{k) is related to the total decoupling temperature by 



(6.7) 



This is because the neutrino interaction rates are proportional to and fcph oc 1/a cx T, hence 

r^/cjph oc T*. 

To compute the damping of acoustic oscillations in the radiation fluid by neutrino diffusion, we fol- 
low Weinberg |95|. For a radiation fluid, shear viscosity is dominant, bulk viscosity vanishes and heat 
conduction is negligible. The shear viscosity is given by 

Wise = ^ X! ''"-'^"i ' '^^■^^ 

i 

e^i denotes the energy density of a neutrino species, t^. is the typical collision time. In the subhorizon 
limit the Navier-Stokes equation and the continuity equation give 

5" + ^77viscfc5' + cj2<5 = 0, (6.9) 

Ctot 

a damped oscillator The damping factor for the mode fc at a given conformal time 77 is 



D{k, rj) = exp 



1 

'2 Jo 



(fcph/etot)?yviscfcrf?/ 



(6.10) 



Here the upper limit of the integral is the conformal time r^max = min[77, ?7dcc(fc)], because collisional 
damping of the 7-^^-hadron fluid by neutrinos ceases at the decoupling of the mode k considered. The 
damping per oscillation is largest for ujt^ = cj/F^ « 1, i.e. immediately before neutrino decoupling for 
a given mode. But note that subhorizon modes get strongly damped long before the mode decouples from 
neutrinos, because a weak damping per oscillation is compensated by many oscillations per Hubble time. 

For a first application, we ask what modes are already damped by the time of the QCD transition. At the 
QCD transition, T = T*, the interaction rates for electron and muon neutrinos (resp. antineutrinos) with 
the leptons are the same, F^_. = Fp^ = F^ = Fp^ = 3.1GpT^ [129], since electrons and muons are still 
relativistic. The r-neutrinos interact only via neutral currents with the leptons and have a lower interaction 
rate, T^^ = Fp^ = 0.6G|T^. Evaluating the damping integral, eq. ( 16.1 Oi l, at T = we find that the 
damping factor D{k,Ti,) is < 1/e for {k-p\^/H)T^ > 10"*, which corresponds to A/cdm < 10^^" A/©, i.e. 
acoustic oscillations on these small scales are wiped out before the QCD transition. Therefore no peaks 
in the radiation or in the CDM transfer function can develop below this scale. This small-scale cut-off is 
independent of the bubble separation scale. 

Next we consider T < T^'^'^{H) w 1.4 MeV and evaluate the final damping factor D{k). At the time 
of neutrino decoupling we take a purely radiation-dominated Universe, consisting of only 7, and // to 
evaluate the damping. Muons have disappeared since T and since <Si tn. The interaction rate 
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for electrons is still given by charged and neutral currents, T^^ = Fp^ = l.SGpT^, whereas muon and tau 
neutrinos have the same lower (neutral current) rate, = Fp = Ti^^ = Fp^ = 0.3GpT^ II129I . The 
final damping of a certain scale (expressed in invariant CDM mass) from eq. (I6.10> is 



with the neutrino damping scale A/,,_dmp ~ 1-9 x 10"'' Mq in CDM mass. This is 3 x lO^^A/^'^™ at 
T — lA MeV and corresponds to length scales A = l/30i?H- 

We conclude that any small-scale density (temperature) fluctuations of the radiation fluid are damped 
away by the time of BBN. The large peaks in the radiation fluid that could be generated by the QCD 
transition do not give rise to inhomogeneous BBN. 



7 Isothermal baryon fluctuations 

The nucleation of bubbles in a first-order QCD transition may strongly affect the distribution of baryons 
in the Universe |20 21 1. Isothermal baryon number fluctuations are created at the end of the QCD phase 
transition. These fluctuations are smoothened by baryon diffusion after the QCD transition. As the tem- 
perature of the Universe drops below 1 MeV, neutrons and protons decouple chemically (freeze out). The 
mean free path of the neutrons exceeds that of the protons and protons and neutrons are thus segregated. 
In other words, regions of high baryon density become proton-enriched and regions of low baryon density 
become neutron-enriched. This affects the local nucleosynthesis yields, giving rise to inhomogeneous big 
bang nucleosynthesis (IBBN) 1311 (in contrast to homogeneous BBN L35J). which is sketched in section 



Isothermal fluctuations of baryon number are an inevitable consequence of any first-order QCD transition 
in the early Universe. After the tiny period of supercooling, bubble nucleation, and reheating, the HG 
bubbles grow adiabatically and merge (percolate). For tiny supercooling the percolation happens very 
close to equilibrium, because this process is extremely slow with respect to the QCD time scale. Just 
before the transition is completed, the Universe is filled with a phase of HG in which shrinking drops 
of QGP are embedded. Usually these drops are assumed to be spheres, which is the case if the surface 
tension is reasonably large. In the homogeneous nucleation scenario, the typical distance between QGP 
drops is given by the typical bubble nucleation distance d^uc- The bulk properties of growing HG bubbles 
and shrinking QGP drops (except for the very last moments before the drops evaporate) are governed 
by the hydrodynamic equations for the radiation fluid only, i.e. the baryons do not influence the bulk 
properties. The hydrodynamics of growing HG bubbles and decaying QGP drops has been studied in 
detail in Refs. 1231 1241 [T30l . Without baryons, no relics would remain from the bubble scale today, so let 
me concentrate on the role of baryons during a first-order QCD transition. 

Let me assume for a while that chemical and thermal equilibrium are maintained from the formation 
of the QGP drops until the completion of the phase transition (evaporation of the last QGP drops). In 
chemical equilibrium the baryon number prefers to reside in the QGP, since the deconfined quarks are 
almost massless. The ratio k of the baryon number density in the two phases is a constant for chemical 




(6.11) 



7. 1 Baryon inhomogeneities 
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equilibrium, given by 



3 



9 .2/ KT.) 

w 16 for = 160 MeV , 

evaluated for the bag model with Nf — 3. The mean baryon number density is given by 

ne = (1 - f)nf^ + fnf" , (7.3) 

where / is the volume fraction filled by the HG. It rises from to 1 during the transition. After reheating, 
/ ^ lO"'^. At percolation / ~ 0.7 (for spherical bubbles with a single radius) 1 131 1. The baryon number 
= tibV is conserved in the comoving volume V. From eqs. ( 17. 1> and il.3\ one finds that both 
rig'^^/riB and ng^/yiB grow during the transition: 



QGP _ '^B 

HQ _ 



(7.5) 



The maximum values are reached for / sa 1, where ng^^^/ne ~ k = 0(10) andng*^/nB ~ 1. Of course, 
if equilibrium were maintained throughout the transition, no baryon inhomogeneity would arise from the 
preference of baryons to reside in the QGP. If baryons drop out of equilibrium before the transition is 
completed, baryon inhomogeneities are necessarily produced. 

There are three mechanisms that drive the baryons out of chemical equilibrium: 

• The baryon flux through the surface of the QGP drop might be suppressed. This was elaborated in 
Refs. 1,21. 26. .81,] . Any suppression of the baryon flux results in a deviation from chemical equilib- 
rium. Since the net flow of baryons is directed from inside the QGP drops to the HG, flux suppression 
increases the value of k and baryon number fluctuations are generated. 

However, we do not know the physics of the phase interface so far It is thus unclear how large this 
flux suppression might be. The simplest model is based on phase-space arguments |26|. A more 
sophisticated model attaches chromoelectric flux tubes to the surface at the points where a quark 
penetrates the interface. These flux tubes then break up and produce mainly pions, and rarely nucleons 
or heavier baryons |81 1. Both models give a large suppression of baryon number flux. However, in 
both models it is assumed that no bound states or cluster states exist in the QGP near the critical 
temperature. If, for instance, a significant fraction of quarks in the QGP would form diquarks near the 
critical temperature, the arguments for flux suppression would break down. 

Our ignorance about the physics of the phase interface has been encoded in the so-called filter factor F 
p6|, which is any number between and 1 and gives the probability that a baryon penetrates the 
interface. The flux suppression is strongest and most simple to analyse when F — 0, which means 
that the baryon number resides in the QGP phase. In such a situation the only relevant quantity is 
the typical nucleation distance. Any QGP drop contains a baryon number of about N-q — UBdf-^^^, 
which gives ~ 10^" (10^^) baryons per QGP drop for dnuc ^ 10"^ (10"^)i?H after the annihilation 
of all antibaryons in the QGP drop. The baryon density in the QGP drops reaches nuclear matter 
density (~ 200 MeV/fm'^) at a drop radius of 10 /im (1 mm). This demonstrates that the produced 
inhomogeneities may be enormous. On the other hand it is clear that the filter factor cannot vanish 
exactly, because all baryons would be trapped in the QGP drops (quark nuggets) by the end of the 
QCD transition and no nucleosynthesis could have taken place. 
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• Even if the baryon flux through the phase interface is large enough to maintain chemical equilibrium 
locally, we must know the diffusion length of baryons during the evaporation of the QGP. Only if 
baryons diffuse over scales dnuc can chemical equilibrium be maintained globally. According to 
eq. il.Si . ng*^ « «b/k ^ '^b at the early stages of the transition, and thus the regions where HG 
bubbles nucleate are underdense in baryon number. 

The typical collision time for baryons is Tb ^ l/T^ (collisions with vrs dominate at T^,). The duration 
of the transition is ^ O.ltn- Thus, baryons make N ^ 10^^ collisions during the phase transition. 
Their typical velocity is Wb ~ (2*/'7^p)^^^ and thus they diffuse a distance of Ab-diff' ~ V^WbTb ~ 
10~^^i?H, which is much smaller than dnuc- Thus baryon diffusion alone is not able to keep the mixed 
phase in global chemical equilibrium. 

The consequences of baryon diffusion have been studied by Kurki-Suonio 11321 and have been in- 
corporated into numerical simulations 11331 . Since i?diop ^ Ab-diff initially, chemical equilibrium 
cannot be established globally (not even within a single phase). At the QGP side of the interface, 
the baryon number is piled up in a thin surface layer of dimension rg ~ A^_jj;ff /i?drop 1 132 1. Since 
locally = Kn^'^, more baryons leak into the HG than without diffusion. The baryon density 

grows when the QGP drops shrink, and thus the baryon density in the HG close to the QGP drop 
grows as well. Since the diffusion length is small the baryons that penetrate the interface cannot be 
redistributed in the whole volume. In this way core-concentrated clouds of baryons are formed. 

• The third mechanism is the extraction of entropy once the mean free path of any relativistic species 
exceeds the size of the QGP drops. This happens first for the neutrinos, A^-mfp ^ 10^^i?H, later 
for the charged leptons. A radiative flux is established, which caiTies away entropy but not baryon 
number. In contrast the hydrodynamic flux carries away entropy as well as baryon number. Thus 
baryons are compressed in the QGP drop (it is the ratio tiq/s that is important). This mechanism has 
been found by Applegate and Hogan |21 1 and was further studied in 11321 126II133I . It turns out that 
entropy extraction is most efficient for large initial drop size (i?diop ^ A^_mfp). 

All three mechanisms lead to the formation of baryon number fluctuations by the end of a first-order QCD 
transition. Because the baryons are tightly coupled to the radiation fluid on scales of ©(dnuc), the baryon 
fluid has the same local temperature as the radiation fluid and thus isothermal fluctuations are generated (in 
contrast to isentropic fluctuations). Let me point out that the second and third mechanisms are independent 
from the details of the QCD phase transition, especially from its order. It might therefore be possible that an 
inhomogeneous (temperature fluctuations from cosmological inflation) crossover could introduce baryon 
inhomogeneities as well. The typical length scale would be Agmooth, as in the case of inhomogeneous 
nucleation. To my knowledge, this issue has not been investigated in any detail yet. 

Let me estimate (for a first order transition) the fluctuation amplitude and the fraction of baryons af- 
fected: the most extreme case takes place for the observationally excluded situation of a vanishing filter 
factor {F = 0). For small values of F we still find extremely high baryon densities and a large fraction of 
affected baryons. Given the large uncertainties on the value of F, it is most interesting to ask what happens 
if F = 1. Since the physics of baryon diffusion and entropy extraction are well understood, this case gives 
a lower limit on the fraction of baryons that is concentrated in clouds. 1 will focus on large nucleation 
distances, 1 take dnuc ^ 10^'*i?H, because, as will be shown below, dnuc > 1 m is necessary in order to 
affect BBN. 

The first substantial departure from chemical equilibrium happens when -Rdrop ^ ^i^-mip ^ 10^^i?H- 
For dnuc ^ 10^'*i?H, the volume fraction in QGP drops at i?drop ^ A^-mfp is /qgp ~ 10^^. Since 
at that time Ab-diff <C i?drop the fraction of baryons that is piled up in a surface layer on the QGP side 
of the interface is negligible. This means that the fraction of baryons that is inside the QGP drops at 
neutrino decoupling is ^ k/qgp, according to eq. ilA\ . and therefore at most ^ IO^^A^b baryons are in 
core-concentrated clouds by the end of the QCD transition. 
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It was shown in Refs. I21II132I that for i?drop ^ Ab-diff the baryon number distribution of such a cloud 
is given by 



where i?cioud ~ Aj^-mfp and a is the ratio of the neutrino radiation flux and the hydrodynamic flux, i.e. 



This ratio depends on the temperature difference between the QGP drop and the surrounding HG. The 
temperature gradient is necessary to establish a hydrodynamic flux. This flux originates from a surface 
layer, whereas the neutrino flux originates from the volume of the drop. When the QGP drop shrinks, the 
ratio a grows because AT has to grow to keep both phases close to Ti,. Within the bag model a ~ 10~^ 
at most. The parameters from quenched lattice QCD give a ^ 10^^. The resulting baryon distribution 
J7.6> is very broad and high baryon densities are reached at the centre of the clouds only. The fraction of 
baryons in the clouds has been estimated to be a/3 1 132 1, which is consistent with the estimate ^ k/qgp- 
In general, increasing the nucleation distance reduces the fraction of baryons in high density regions. 

It was shown numerically in Ref. 1 133 1 that enormous baryon densities are reached in the inner 10^ fm 
of the drop. However, it is hard to have a cosmologically significant fraction of baryons residing in these 
inner regions if there is no baryon flux suppression. Large separations of baryon clouds and high baryon 
fraction within these clouds are mutually exclusive, except for significant baryon flux suppression. 

7.2 Baryon diffusion 

The isothermal baryon number fluctuations formed at the end of a first-order QCD transition are smooth- 
ened by two mechanisms, baryon diffusion 1 1341 1 1351 1321 and neutrino inflation Il21lll29l|l36l . The latter 
is important for baryon clouds with high 77 = ne/n^ > 10^^. The pressure of these baryon clouds is 
in equilibrium with the regions of lower baryon density. Thus the temperature of the baryons could be 
below the mean temperature, because the baryons give an additional contribution to the pressure besides 
the radiation quanta. If the size of the cloud is smaller than the mean free path of neutrinos, heat is 
transported into the baryon clouds. This leads to an increase in pressure inside the cloud and thus to an 
expansion (inflation) of the cloud. Neutrino inflation provides an upper limit 77 < 10~* for clouds with 
^cioud < 1 m at the time of the QCD transition. In contrast to neutrino inflation, baryon diffusion is an 
efficient smoothing mechanism for any value of -q. 

After the annihilation/decay of pions the baryon mean free path increases (below T ^ AQ MeV). At tem- 
peratures above the decoupling of the weak interaction 1 MeV) protons and neutrons are in chemical 
equilibrium and thus nucleons frequently change isospin. A nucleon spends a fraction X^{T) = n^/n-Q = 
1/[1 + exp((5/r)] of its time as a neutron, where Q = 1.29 MeV is the neutron-proton mass differ- 
ence. The mean free path of neutrons is much larger than the mean free path of protons. Thus at T > 1 
MeV the diffusion length of nucleons is given by Aw-diff (0 ^ \/ K-i-mfpVnXnt, where An-mfp is the 
neutron mean free path and Vn ~ {T/m-^Y^^ is the thermal velocity of neutrons. The neutron mean free 
path is determined from neutron-electron scattering (the neutron has an electric dipole moment) and from 
neutron-proton scattering. At 1 MeV the nucleon diffusion length is « 200 m 1 134 , 32 1. This corresponds 
to a comoving distance at the QCD transition of ~ 1 m or ~ 10^'*i?H- 

After weak decoupling (at ^ 1 MeV) neutrons diffuse much larger distances than protons. The proton 
diffusion is determined by electron-proton scattering, which gives a diffusion distance much smaller than 
that of the neutron. The proton diffusion at the time nucleosynthesis starts (at 0.1 MeV) is completely 
negligible compared with the nucleon diffusion before weak decoupling. Thus Ap_diff ^ 3 km at the start 
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ofBBN and 



Ap-diff 10 "^-Rh comoving at T-^. 



(7.8) 



The neutron diffusion length is An-diff ^ 200 km at the start of BBN or 



A, 



•n-diff 



10 ^i?H comoving at T*. 



(7.9) 



The typical scale for baryon inhomogeneities are given by the typical bubble nucleation distance, dnuc. 
Thus all baryon inhomogeneities on scales dnuc < Ap_diff are washed out by baryon diffusion by the 
time BBN starts. Most interesting are baryon inhomogeneities on scales Ap_difi < c'nuc < An-diff. The 
situation dnuc > A^-diff seems to be highly disfavoured, according to the discussion in section |5] Note 
that i?H « 6 X lO'^ km at T = 0.1 MeV, thus Ap_diff 10"'^i?H at the start of BBN. 

The different diffusion lengths lead to the segregation of protons and neutrons. Neutrons escape from 
the high baryon density clouds, while the protons remain in the cloud. Thus the surrounding cloud becomes 
neutron-rich, whereas the baryon cloud itself becomes proton-rich. This provides inhomogeneous initial 
conditions for nucleosynthesis. 

There is no damping mechanism for baryon inhomogeneities on scales larger An-diff. However, it is 
hard to imagine any mechanism that could have generated such fluctuations. As shown above, a first-order 
QCD transition would give rise to inhomogeneity scales just above Ap_diff . Especially, if we assume that 
cosmological inflation took place, baryogenesis must have happend after inflation and thus any baryon in- 
homogeneities must have been produced during the radiation-dominated epoch. Besides the cosmological 
QCD transition, there is no other event before BBN that possibly could have generated isothermal baryon 
fluctuations at large enough scales. If one assumes that the QCD transition does not give rise to any relics, 
the baryon-to-entropy ratio is constant not only in time, but also in space. 



The baryon segregation, which leads to an inhomogeneous neutron-to-proton ratio, and inhomogeneities 
in the entropy per baryon both lead to inhomogeneous nucleosynthesis 131. .138.1 . Let me first discuss what 
happens to the element abundances in a homogeneous patch when we change rj and Xn- 

rj In homogeneous BBN, rj is the single free parameter of the theory. Based on measured deuterium 
abundances, we obtain -q = (5.6 ± 0.5) x 10"^° or cjb = 0.0205 ± 0.0018 EH, consistent wifli the 
value obtained from CMB experiments [36" "3^ . The key role of 77 is to determine the time when the 
nucleosynthesis starts. The first step in nucleosynthesis is the generation of deuterium (p H- n ^ d H- 
7). Since deuterium has a small binding energy {k, 2.2 MeV) it is very easy to photo-dissociate it. The 
condition that deuterium does not see photons with energy above 2.2 MeV is rj-^ exp(2.2 MeV/T) < 
1, which gives the temperature at the beginning of nucleosynthesis. The higher 77, the earlier BBN 
starts, which means that baryons can be burned more efficiently. For higher 77, more ^He and less D 
are produced. Typical mass fractions^ in homogeneous BBN are 0.25, 5 x 10~^, 3 x 10~^, 3 x 10~^ 
for ^He, D, ■^He, and ^Li, respectively. For very high baryon densities, i.e. -q ^ 10~^, heavy elements 
such as ^^C, ^■^C, and ^"'N could be generated with mass abundances of 0(10^^) |4i ]. 

The neutron-baryon fraction determines the mass abundance of ''He, given by = 47i4jjc/«-B ~ 
2Xn. Here I used the approximation that all neutrons eventually end up in ^He. For homogeneous 
BBN X^ w 1/8 and thus Y-p k, 0.25 from this simple argument. An inhomogeneous distribution of 
X„ could change the local ''He abundance. Of course, the other abundances change as well. 



Note, usually for D, ''He and ^Li the number fractions are given. Here I quote mass fractions for all elements. 



7.3 Inhomogeneous nucleosynthesis 



adp header will be provided by the publisher 



39 



In order to produce a global change in the elements abundances, non-Unear fluctuations in 77 and/or are 
needed. Linear fluctuations would yield the same average abundances as homogeneous BBN. In addition 
the IBBN yields depend on the geometry of the baryon fluctuations. 

Originally, the interest in IBBN arose because it was thought that i^B = 1 was possible |134|. It 
was shown that even with IBBN a baryon-dominated, flat Universe is impossible, as it means a severe 
overproduction of ^Li 1 137„ .26„.138,l . Nevertheless there are interesting, probably in the future observable, 
consequences from IBBN: the precision measurements of the primordial abundances will need to take 
effects from IBBN into account 1 33 1 and, more speculative, higher mass elements (C, N, O) might have 
been generated during IBBN at an amount that could be observable 141 J . 

The following conditions are necessary in order to have significant effects from IBBN: 1. dauc > 
Ap-diff, 2. 77high/?7iow should be large, and 3. the fraction of baryons in the high density region has to 
be well above the ^Li abundance. Heavy elements could be produced, leaving the ''He abundance un- 
changed, if this fraction is small and rjhigh is at its upper limit ^ 10^"'. On the other hand, significant 
corrections to the light element yields are obtained if the fraction of baryons in high density regions is 
large. In a recent work, Suh and Mathews |32J found that two fluctuation scales are consistent with the 
observed abundances: < 1 m (homogeneous BBN) and 10-40 m, measured at T*. It seems that the nu- 
cleation distances obtained from heterogeneous and inhomogeneous nucleation are consistent with these 
constraints on IBBN. 

8 Cold dark matter 

The intention of this section is to highlight the differences between various CDM candidates. By definition, 
CDM has negligible pressure and speed of sound when the formation of structure starts (at t ~ tcq). Before 
matter-radiation equality, inhomogeneities in the density of CDM grow logarithmically, whereas thereafter 
their growth is linear in the scale factor. This is the same for all CDM candidates, and thus they cannot 
be distinguished at mass scales above Meq = 1.32 x lO'^^(O.15/a'in)^-M0, corresponding to a comoving 
distance Rcq ~ 91(0.15/win) Mpc today. However, the astrophysics of CDM depends on the nature of 
CDM at scales well below Rcq. Besides behaving like a non-relativistic fluid at large scales, a good CDM 
candidate should be inert to electromagnetic and strong interactions. Its interactions should not be stronger 
than the weak interaction. 

Here I pick three examples: the neutralino, the axion and primordial black holes. The neutralino is 
an example of a WIMP. The general features of a WIMP are that it is massive (m ^ T^q) and interacts 
weakly. The axion interacts weaker than the weak interactions of the standard model and it is very light. 
Nevertheless, it is a CDM candidate, because it oscillates coherently and axion strings decay into non- 
thermal axions. Primordial black holes are very heavy on the particle physics scale and their sole interaction 
is gravity. Before discussing these candidates in more detail, let me mention some of the other candidates 
that have been proposed (strangelets, quark nuggets and QCD-balls are discussed in section|2ji. 

Another CDM candidate is a heavy neutrino (4th generation) with < nii, < 1 TeV 1 139|, where 

only the upper mass bound gives cosmologically relevant CDM. Instead of the neutralino, the axino, the 
supersymmetric partner of the axion, could be the lightest supersymmetric particle. Axinos (m^ < 300 
eV) would be a non-thermal relic, generated by the decay of neutralinos at T ^ 10 MeV 11401 . Other 
non-thermal relics are WIMPzillas |141|, particles that are so heavy that they can never be in thermal 
equilibrium. They would be generated at the end of cosmological inflation. 

8. 1 Neutralinos and other WIMPs 

Let me start the discussion with the lightest supersymmetric particle [1421, which is one of the most popular 
CDM candidates nowadays. In the minimal supersymmetric standard model (MSSM) |143| the lightest 
supersymmetric particle most Ukely is the neutralino (when it is assumed to be stable). The neutralino is 
a mixture of the bino, neutral wino, and the two neutral higgsinos. The lightest one is expected to closely 
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Fig. 8.1 Feynman diagrams contributing to elastic neutralino-lepton scattering; x is the neutralino, I is a lepton, and 
Il,r denote the corresponding left- and right-handed sleptons. 



resemble a bino [144]. Assuming high slepton masses (> 200 GeV) and the unification of gaugino mass 
parameters Mi and AI2 at the GUT scale, allows to derive constraints on the mass of the lightest neutralino 
from the LEP 2 experiments. The result is to^ > 39 GeV |5 1. Adding constraints from cosmology (CMB 
fits to a ACDM model) and from 6 ^ 57 decay, suggests that 100 GeV < < 500 GeV'° |145|. 

It is essential to distinguish between the chemical freeze-out and the kinetic decoupling of neutralinos: 
the chemical freeze-out determines the amount of neutraUnos today and occurs when the annihilation 
rate of neutralinos drops below the Hubble rate, Tann/H < 1. Soon after the neutralinos become non- 
relativistic, the rate for neutralino annihilation, Fann = (wcrann)f^X' suppressed by the Boltzmann factor 
in the number density of the neutralinos, n-j^ ~ (rn^T)^/^ exp(— m-^^j/T). The freeze-out temperature 
of the neutralino fT42l is Ted ~ m^/20 > 5 GeV, and neutralinos are chemically decoupled at the QCD 
transition. We do not discuss the physics of neuralino freeze-out in any detail here. Besides the annihilation 
of neutralinos, it is important to take into account coannihilation channels (e.g. a neutralino and a sfermion 
could 'annihilate' into ordinary particles). If the masses of the lightest neutralino and other sparticles are 
similar, coannihilation plays an important role. Also contributions from poles and thresholds have to be 
treated with care. For details we refer the reader to refs. 1145111461 . 

Kinetic decoupling, in contrast, is determined by the elastic scattering between neutralinos and the 
dominant radiation fluid. The interaction rate for elastic scattering is Fd — r^^Jj = {vaci)n, where n ^ 
is the number density of relativistic particles, e.g. electrons or neutrinos. The relevant diagrams for elastic 
neutralino lepton scattering are shown in figure lO] Other contributions have been demonstrated to be 
of subleading order [1491 . If the neutralino is a pure gaugino, there is no contribution from Z exchange. 
An order of magnitude estimate shows that {ac\v) ~ GpT^ I147II13I . similar to elastic neutron-neutrino 
scattering [1481. 

One has to distinguish between the regime of perfect kinetic coupling, i.e. neutralinos tightly coupled 
to the radiation fluid, an intermediate regime where the neutralinos scatter elastically but the number of 
collisions is not sufficient to drag them along dissipationless, and the regime of free streaming, Vc\/ H < 1, 
which is roughly at T < 1 MeV, since the neutralino interacts weakly. 
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The upper limit holds for tan /3 < 30(45) for fi < 0(> 0) in the parameter space of the constrained MSSM. 
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Fig. 8.2 Kinetic decoupling as a function of the universal sfermion (slepton) mass for bino masses £ 
{50, 100, 200} GeV (from bottom to top) 1,131 . 

Let us estimate the regime where the neutralinos belong to the perfect radiation fluid. Perfectness of 
a (dissipationless) fluid refers to an external time scale, here uj^^ of an acoustic oscillation or the Hubble 
time. A fluid behaves as a perfect fluid if ujT^ciax ^ 1, i-S- the external time scale is larger than the 
relevant relaxation time, the fluid is continually in local thermal equilibrium, and no entropy production 
occurs. For the coupling of neutralinos to the radiation fluid the relaxation time is given by Trciax = NtcoW, 
where N is the number of collisions needed to completely change the momentum of the neutralino due to 
collisions with the radiation fluid. The momentum transfer at a collision of a lepton with the neutralino 
is of order p\ ^ T. The kinetic energy of the neutralino is given by equipartition, p^/2mj^ ~ T, hence 
its momentum is ~ ^ra^JT . The fractional change of the neutralino momentum from one collision at 
the QCD transition is S-p^j-p^ ^ ^jT/m^ <C 1. After N collisions the total rms change of momentum is 
[5Px/Px)rn\s ^ v^A^T/mj^. Local thermal equilibrium is obtained if the bulk motion of the neutralinos is 
governed by the leptons, i.e. the fractional change of the neutralino momentum is of order 1. The number 
of collisions needed to completely change direction is ~ ™x/^ ~ ^^^^ for ra^ — hQ GeV and T = Ti,. 
The colhsion time is given by the weak interactions rate ~ t^^Ji ^ GpT^, and the relaxation time is 
given by 

T-reiax - Nt,ou ^ VV X IQ-'in, (8.1) 

with N ~ m^/T. If the relaxation time is compared with the frequency of the acoustic oscillations, one 
finds that the condition for a perfect fluid, Ci^Ticiax <C 1, is satisfied, at the QCD transition, for scales A > 
^x-dcc{T-t) ~ 10^^i?H (nij^ = 50 GeV). Hence the neutralinos on these scales are part of the radiation 
fluid at the QCD transition. Below this scale, the neutralinos cannot follow the acoustic oscillations. On 
the other hand on the Hubble scale the perfect kinetic coupling of neutralinos to the radiation fluid stops 
when the required relaxation time becomes more than a Hubble time, Trdax > ^h- This gives a temperature 
of T^-doc ^ 10 MeV The kinetic decoupling of neutralinos has been studied in detail in fl31. In figure 
18.21 the dependence of the kinetic decoupling temperature on the neutralino mass (here a bino) and the 
slepton mass is shown. Down to the decoupUng temperature neutralinos on the Hubble scale belong to the 
radiation fluid. 
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The kinetic decoupling of other WIMPs, such as heavy neutrinos, also happens way below T*. Therefore 
the neutralino or a heavy neutrino would be tightly coupled to the radiation fluid during the QCD transition. 



The axion has been introduced to solve the strong CP problem and it is a popular CDM candidate 11501 . 
Direct searches in the laboratory, constraints from astrophysics and from cosmology allow the axion 
mass to be 10^^ - 3 x 10"'^ eV". Axions could be the dominant matter today if their mass is small, 
i.e. TOa ^ 10^^ eV, which corresponds to a breaking of the Peccei-Quinn (PQ) symmetry at the scale 
/pQ ~ 10^^ GeV. These axions could be produced coherently by an initial misalignment of the axion field 
and by the decay of axionic strings (the latter is only possible if the reheating temperature after inflation 
lies above the PQ scale). 

The initially misaligned axion field starts to oscillate coherently when the axion mass has grown to 
rna(Ti) ~ 3iJ(Ti), where Ti ~ 1 GeV LI 52.1 . Below Ti - 1 GeV the oscillating axion field evolves 
as CDM. This can be seen as follows: the energy density of a scalar field is p = (^^/2 + V{(p) and the 
pressure is p — "yi^/S — V{(p). For time scales large with respect to the period of the coherent oscillations 
Ta ^ TTia ^ ~ lO^^^H, the virial theorem can be applied to give (pa) — 0. Thus, a coherently oscillating 
scalar field, like the axion field below Ti, can be modelled by a pressureless fluid. 

The axion field is inhomogeneous on scales larger than Rb_{Ti), if, as is usually assumed, the reheating 
temperature after inflation lies above /pq. This leads to inhomogeneities with 5pa/ Pa. = in the axion 
density. These inhomogeneities evolve into axion miniclusters of typical mass M^c ^ lO^^^M© and 
typical radius i?inc ~ 10^ m ^ QARq ^ au [1531. If the reheating temperature was below the PQ 
scale, the axion field at Ti was homogeneous on superhorizon scales, apart from small, quantum-induced 
fluctuations of the axion field during inflation. 

To summarize, the axion is a CDM candidate that is kinetically decoupled from the radiation fluid at the 
QCD scale. This is also the case for WIMPzillas and for primordial black holes. 

8 . 3 Primordial black holes 

A further possibility for CDM that decoupled kinetically long before the QCD transition is primordial black 
holes (PBHs) [128. 126.1 produced before the QCD transition and therefore with masses A/bh ^ IMq. 
In order to survive until today, PBHs should have A/bh > lO^^g « 10^ ^^Af©. PBHs in the range from 
\Q~^^Mq to lO~^^Af0 would radiate too strongly to be compatible with 7-ray observations [1 05111541 . 
The production of PBHs arises naturally from hybrid inflationary scenarios 1 155. .154, .127..156J . 

When these black holes are formed, they are comoving with the radiation fluid, but their velocity is 
redshifted and soon becomes completely negligible. PBHs should not overdose the Universe, so that just 
a tiny fraction of energy in a comoving volume is allowed to collapse to black holes. This fraction is given 
by / ~ apBH/flcq , where opbh is the scale factor at PBH formation and Ooq at matter-radiation equality. 
Their density at the QCD epoch is 



In every Hubble volume PBHs make up a pressureless gas if ripBH(7*) S> ^^^^ ^pbh > lO^Tt ~ 

1 TeV This corresponds to PBHs with mpsn < lO^^Af©. PBHs smaller than this mass make a kinetically 
decoupled fluid at the QCD horizon scale. 

^' A stronger limit is obtained in fisn . i.e. nia, > 3.7 X 10~* eV. 



8.2 Axions 




(8.2) 



at most, where Tpbh is the temperature at which PBHs form with mass 




(8.3) 
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9 Density fluctuations in CDM 

9 . 1 Kinetically decoupled CDM 

CDM is assumed to be a major component of the Universe today. At the time of the QCD transition, 
however, the contribution of CDM to the total energy density was negligible: ecdm/eiad = a /acq ~ 
10 Here, we consider a type of CDM which is non-relativistic {p <C e) at the QCD transition and 
which is only coupled via gravity to the radiation fluid. Examples are axions or PBHs. CDM moves in 
the external gravitational potential provided by the dominant radiation fluid. During a first-order QCD 
transition, the big amplifications of the density contrast in the radiation fluid (Jiad (see section lSTl leads to 
a big amplification in the gravitational potential. The CDM is accelerated to higher velocities at the end of 
the transition. Therefore, we also get peaks and dips in the cold dark matter fluid, as has been discussed in 

The initial conditions for CDM are obtained assuming adiabatic perturbations, i.e. the entropy per cold 
particle is unperturbed (5 (srad/?tcdm) = 0. Since ecdm oc and trad oc 1/a'*, the adiabatic initial 

conditions for (5cdm can be written Jcdm = (3/4)(5rad- The initial fluid velocities are equal. 

The subhorizon evolution of CDM in a purely radiation-dominated Universe is just inertial motion, as 
can be seen from the Euler equation. At leading order in a; = k/H, the gravitational force can be neglected 
and the subhorizon evolution of the CDM velocity is obtained: 

4dm = C-. (9.1) 

X 

The velocity of CDM in a radiation-dominated Universe just redshifts to zero on subhorizon scales; C is an 
integration constant of order A-m- The corresponding evolution of the density contrast Scdm follows from 
the continuity equation, 

Scdm^Clnx + D. (9.2) 

This logarithmic growth of Scdm [1571 can be seen in fig. l9.1l before and after the transition. The shape of 
(5cdm can also be seen in the transfer functions of figs. l6.1l and l6.2l on scales above the horizon scale Mi. 
The major amplification effect comes from a higher velocity at the end of the transition, which leads to 

Scdmiri) = i'f"'kv+ In (^-jL^ + 5f''\ (9.3) 
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The amplification of the density contrast during the transition is neghgible with respect to the enhancement 
of the velocity. This velocity enhancement during the transition leads to an additional logarithmic growth 
of ^cdm after the transition: 



krjj^ In 



rcdm 



(9.4) 



This additional logarithmic growth of (5cdm is most clearly seen in fig. 19. II 

Let us compare the 'QCD peaks' in the CDM spectrum with the CDM spectrum without phase transi- 
tion. In the limit ecdm/ciad <C 1, without transition, we find (x — k/Ti) in the subhorizon limit: 



^cdm — 




(9.5) 



which fixes the constants C and D in eqs. ( 19.11 1 and (|^}. 
We are now able to calculate the enhancement factor 



6^ 



transition 



Jno transition 



(9.6) 



cdm 



The enhancement at the smallest scales Mcdn 
intermediate scale Mcdm ~ 10~^^Mq reads 



10 Mq to which our calculation applies and at some 



600(20) bag model 
10(2) lattice fit 



forl0"^°(10"^^)Mo 



(9.7) 



This shows that both models lead to large enhancements of the CDM density contrast at equality for small 
enough scales. For the lattice fit, the enhancement becomes important at scales below Afcdm ^ 10~^^Mq. 

9.2 CDM clumps 

CDM in the form of axions or PBHs is not subject to damping as the radiation fluid, thus the peaks in CDM 
will survive until structure formation starts. The free streaming scale of CDM is way below our smallest 
scales, because the initial velocity of axions or PBHs is completely negligible. An overdensity of CDM in 
the form of axions or PBH (or any other matter that is kinetically decoupled at the QCD epoch) decouples 
from the cosmic expansion when its density contrast becomes non-linear, {5p/ p)r ~ 1 (condition for turn 
around). It collapses and virializes by violent gravitational relaxations and forms a clump of CDM. For 
spherical collapse the final viral radius is half of the radius at turn around 1 158 1. In the following, I present 
an updated version of the discussion in 1 52 1. 

We take a COBE |93| normalized spectrum and allow for a tilt |n — 1| < 0.2. During the radiation- 
dominated regime, {Sp/p)k for CDM continues to grow logarithmically. In {5p/ p)\ another logarithm 
comes from summing modes up to k ^ 1/R, where R is the radius of the window function. The en- 
hancement factor £ of CDM density fluctuations during the QCD transition has been obtained in eq. ( 19. 6> . 
Putting these factors together we obtain the amplitude of CDM perturbations of size R at the time of 
equality: 



f \ cdm 

P Jr. 



2 X 10-* I — 

fco 



£k{tcq) 



(9.8) 



where ko denotes the comoving wave number of the mode crossing the horizon today. In the following we 
assume uj^ = 0.15, and thus z^q w 3600. For a CDM mass of lO^^^ (lO^^o^^j^ ^ jjjj of n - 1 = 0.2 
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gives a factor « 16 (23). The logarithms contribute another factor « 98 (144). The enhancement factor 
has been calculated in section lgTl to be £{tcq) ~ 2(10) for the lattice QCD fit. Looking at Scr (= 3 standard 
deviations) peaks, we find without tilt (n — 1 = 0): 

c \ 3(T,n— 1 

— ) (icq) « 0.1(0.9) forAfeiump- 1O-1'(1O-2O)M0 . (9.9) 

P / R 

This implies that these clumps go non-linear at Zn\ ~ 360 (3200) and collapse to clumps of radius 
i?ciump ~ 23(O.O5)i?0. With tilt (n - 1 = 0.2), we find 

c \ 3(T,n— 1.2 

— ) (icq) « 2(20) forMciump = 1O-1'(1O-2")M0 . (9.10) 

P / R^ 

These over-densities start to collapse even before tcq, because CDM locally starts to dominate at ~ 
2(20)req. This leads to clump radii of i?ciump « 1.2(0.004)i?Q. 

We conclude that the peaks in the CDM spectrum lead to clumps of masses 10^^° - 10~^"AfQ. To- 
day, these clumps would have a density contrast of 10^° - 10^^, where the lower value corresponds to a 
IO^^'^A/q clump from an untilted CDM spectrum, the larger value is for a lO~^°Af0 clump from a tilted 
CDM spectrum. The evolution of these clumps in the late stages of structure formation remains to be 
investigated (disruption, mergers, etc.). 

For larger enhancement, e.g. if it should turn out that the latent heat is larger than the value from present 
lattice QCD calculations, more compact clumps are possible. These could be subject to femto-lensing 
|97|. With the values of the lattice fit l!5Tl l52l . the CDM clumps are not compact enough to lie within the 
Einstein radius, which is i?E ~ 0.02i?Q for a 10~^^Mq clump. 

The clumping of CDM changes the expected reaction rates for some dark matter searches, because some 
of the rates depend on the space-time position of the detector, star, or planet. Thus experiments looking for 
axion decay in strong magnetic fields |96 150 1 actually yield upper limits on the local axion interaction 
rate. It seems these experiments tell us that the Earth is not sitting in an axion cloud, if such clouds existed. 



9 . 3 Kinetically coupled CDM 

Let me now turn to CDM candidates that belong to the radiation fluid at some point and decouple kinetically 
long after freeze-out. The most important example here is the neutralino, but the physical processes are the 
same for any WIMP that has been in thermal contact with the radiation fluid at some point. I follow the 
discussion of fl3\ . which are in agreement with the recent findings in II6II . 

Collisional damping 

During the process of kinetic decoupling the neutralinos acquire a finite mean free path. Density in- 
homogeneities on scales of the diffusion length are damped by the mechanism of collisional damping. It 
is convenient to describe the CDM as an imperfect fluid. We have shown in [131 that the coefficient of 
heat conduction vanishes at the leading order in T/m^. Thus the dominant contribution to collisional 
damping comes from bulk and shear viscosity. Since the energy of the CDM fluid can be transferred to the 
radiation fluid, which acts here like an inner degree of freedom for the CDM particles, the bulk viscosity 
does not vanish. Nevertheless, the radiation fluid can be treated as a perfect fluid since eiad ^ Ecdm at 
kinetic decoupling of the neutralinos. We calculated the relevant coefficients of transport from the kinetic 
theory in 1131 . At linear order in the relaxation time the coefficients of shear and bulk viscosity become 
^visc ~ ^ x^'Treiax and Cvisc ~ (5/3)77visc, respectively. 

The density inhomogeneities in CDM are damped exponentially below the scale Afd due to viscosity 

113 ca 
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Fig. 9.2 Acoustic damping scale as a function of the bino mass for sfermion masses Alp £ {150, 200, 300, 400} 
GeV (from bottom to top) 1 13 



In figure l9^ we plot the damping mass M^_dmp as a function of the neutralino mass for various values of 
the slepton mass. The damping ( 19.111 1 provides a small-scale cut-off in the primordial spectrum of density 
perturbations in neutralino CDM. 

Free streaming 

Once the temperature in the Universe drops below Tkd the rate of elastic scatterings is not high enough 
to keep the neutralinos in thermal equilibrium with the radiation fluid. The neutralinos enter the regime of 
free streaming. This process continues to smear out inhomogeneities, since the individual neutralinos do 
not move coherently. From kinetic theory one can show that the damping due to free streaming goes as 



Px 



oc exp 



kd 



2m^ 



If 



InM — 

akd 



exp 



AJfs(Q) 
M 



2/3' 



(9.12) 



The mass scale of damping from free streaming Mfs is written as a function of the cosmic scale factor a. 
In the radiation-dominated Universe, the damping scale grows logarithmically with the scale factor. This 
calculation agrees with the estimate of the free streaming scale from the free streaming length as presented 
in itOl up to a numerical factor (27r/ ^6)-^ w 17. (We previously underestimated the free streaming mass 
by that factor) The ratio 



X-dmp 



3 flkd 



H 3 



(9.13) 



exceeds unity for a > 2.20^^. Free streaming thus starts to dominate the damping from collisional damping 
once the Universe has doubled its size after kinetic decoupling. It is interesting to evaluate ( 19.131 1 at the 
time of matter-radiation equality, since this is the moment when CDM density perturbations start to grow 
linearly with the expansion. For a kinetic decoupling temperature Tkd — 40 MeV and for a;,n = 0.15 we 
find Mfs(aoq)/Mx-dmp ~ 1-3 x 10'', thus Mfe(acq) w 8 x IQ-^Mq for = 150 GeV. 

Typically the free streaming mass at the time of equality is of the order of IO^^-IQ^^Mq, which is 
in striking contrast to claims in the literature (see e.g. 1 160|) that the minimal mass for the first objects 
would be ^ 10^^'^(150 GeV/TO^)'^AfQ. The huge difference with our result comes mainly from the false 
assumption that kinetic decoupling occurs simultaneously with chemical decoupling. 
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Fig. 10.1 History of the first second of the Universe and the epoch of BBN. A sequence of thermodynamic transitions 
takes place: the electroweak transition (~ 10 ps) is followed by the QCD transition (~ 10 /xs) and by annihilation 
('^ 100 s). With respect to particles of the standard model, several important processes occur: baryons and antibaryons 
annihilate (~ 1 ms), neutrino oscillations set in at ~ 0.1 s, neutrinos decouple and at the same time the neutron-to- 
proton ratio freezes out (~ 1 s). During the epoch of BBN (~ 3 min) this ratio changes slightly due to neutron decay 
(r — 886 s). Regarding the cold dark matter candidates, WIMPs freeze out at ~ 10 ns and decouple kinetically at ~ 1 
ms. The axion mass is switched on close to the QCD transition at ~ 1 /xs. 



To summarize, collisional damping and free streaming smear out any primordial density inhomo- 
geneities in neutralino CDM below ^ 10^^ Af©. This implies that there is a peak (subhorizon CDM 
density perturbations grow logarithmically during the radiation epoch) in the power spectrum close to the 
cut-off and therefore we have found the minimal mass for the very first objects, if CDM is made of neu- 
tralinos. This result does not depend in a strong way on the parameters of the supersymmetric standard 
model. Looking at Na overdense neutralino regions, we find that they go nonlinear at z^i ~ 36 for a flat 
spectrum of primordial fluctuations. 

According to the picture of the hierarchical formation of structures, these very small first objects are 
supposed to merge and to form larger objects, eventually galaxies and larger structures. It is unclear 
whether some of the very first objects have a chance to survive. CDM simulations show structures on all 
scales, down to the resolution of the simulation fl62"163V However, the dynamic range is not sufficient to 
deal with the first CDM objects, so the fate of the first CDM clouds is an open issue. A cloudy distribution 
of neutralino CDM in the galaxy would have important implications for direct and indirect searches for dark 
matter. A first attempt to study these consequences in detail has been made recently 1 161 1. It is estimated 
that a small fraction of these neutralino clouds survives tidal disruption and that these are probably enough 
to dominate a neutraUno annihilation signal in our galaxy. 



10 Summary and concluding remarks 

The general conclusion from this work is that many important cosmological issues are linked to the first 
second of the Universe, especially the nature of dark matter and the formation of nucleons. In figure [TOJl 
a summary of the history of that epoch is displayed. 

It might be that no relics from that epoch survive and it will therefore be hard or even impossible to 
probe the first second directly. Nevertheless, it is mandatory to study it in detail, since this is essential to 
answer the following questions: (i) Why are the initial conditions for the BBN (from observations) very 
close to homogeneous? (ii) What is the small-scale structure of CDM and what is the implication of the 
CDM small-scale structure for dark matter searches? (iii) Can we exclude that CDM is baryonic? (Are 
quark nuggets or other QCD-relics excluded)? 

Let me finally stress that this review represents my personal point of view and, there are certainly more 
issues than discussed here that are relevant for a complete understanding of the early Universe. 
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a 


scale factor 




B 


depending on context, baryon number or bag constant 


Cs 


speed of sound 




d 


distance 




f 


depending on context, frequency or volume fraction 


g 


effective number of relativistic helicity degrees of freedom 


H 


Hubble rate 




h 


depending on context, = Hq/{1QQ km/s/Mpc) or amplitude of gravitational waves 


k 


comoving wave number 




L 
I 


lepton number 
latent heat density 




M 


masses of macroscopic objects (with some noted exceptions) 


m 


particle masses 




Nf 


number of quark flavours 




n 


number density 




P 


pressure 




R 


radius 




s 


entropy density 




T 


temperature 






transition temperature 




t 


cosmic time 




V 


velocity 




z 


redshift 




r 


interaction rate 




e 


energy density 






bulk viscosity 




n 


depending on context, = ns/ 


or conformal time 


71 • 

'/Vise 


shear viscosity 




A 


wave length 






chemical potential 







mass density 




(J 


depending on context, surface tension or cross section 


T 


time scale 




f2 


fractional energy density 






= nh"^ 




Table 1 


List of symbols. 







today eq 


matter-radiation equality 


ann 


annihilation fs 


free streaming 


b 


baryons gw 


gravitational waves 


c 


critical H 


typical scale set by Hubble expansion 


cd 


chemical decoupUng kd 


kinetic decoupUng 


cdm 


cold dark matter m 


matter = baryons + CDM 


coll 


collision mfp 


mean free path 


dec 


decoupling nuc 


nucleation 


defl 


deflagration ph 


physical 


diff 


diffusion rad 


radiation 


dmp 


damping rms 


root mean square 


el 


elastic sc 


supercooUng 



Table 2 Meaning of suffixes. 



